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ABSTRACT 

We calculate the structure and neutron content of neutrino-heated magnetohydrodynamic winds driven from 
the surface of newly-formed magnetars ("proto-magnetars") and from the midplane of hyper- accreting disks, 
two of the possible central engines for gamma-ray bursts (GRBs) and hyper-energetic supernovae (SNe). Both 
the surface of proto-magnetars and the midplane of neutrino-cooled accretion flows (NDAFs) are electron de- 
generate and neutron-rich (neutron-to-proton ratio n/p ^> 1). If this substantial free neutron excess is preserved 
to large radii in ultra-relativistic outflows, several important observational consequences may result. Weak in- 
teraction processes, however, can drive n/p to ~ 1 in the nondegenerate regions that obtain just above the 
surfaces of NDAFs and proto-magnetars. Our calculations show that mildly relativistic (Lorentz factor V < 10) 
neutron-rich outflows from NDAFs are possible in the presence of a strong poloidal magnetic field. How- 
ever, neutron-rich winds possess a minimum mass-loss rate that likely precludes simultaneously neutron-rich 
and ultra-relativistic (r > 100) NDAF winds accompanying a substantial accretion power. In contrast, proto- 
magnetars are capable of producing neutron-rich long-duration GRB outflows ~ 10-30 seconds following 
core bounce for sub-millisecond rotation periods; such outflows would, however, accompany only extremely 
energetic events, in which the GRB + SN energy budget exceeds ~4x 10 52 ergs. Neutron-rich highly relativis- 
tic outflows may also be produced during some short -duration GRBs by geometrically thick accretion disks 
formed from compact object mergers. The implications for r-process nucleosynthesis, optical transients due to 
non-relativistic neutron-rich winds, and Nickel production in proto-magnetar and NDAF winds are also briefly 
discussed. 

Subject headings: stars: neutron — stars: winds, outflows — supernovae: general — gamma rays: bursts — 
stars: magnetic fields — accretion, accretion disks — nuclear reactions, nucleosynthesis, 
abundances 



1. INTRODUCTION 

The rapid variability and large energies that characterize 
cosmological gamma-ray bursts (GRBs) strongly implicate 
stellar-mass compact objects as their central engines. In- 
deed, the association of several long-duration GRBs (LGRBs) 
with Type Ibc supernovae (SNe) suggests that LGRBs result 
from relativistic outflow accompanying rapid accretion onto 
a newly-formed black hole (a "collapsar"; Woosley 1993; 
MacFadyen & Woosley 1999) or the spin-down of a newly- 
formed magnetar (e.g., Usov 1992; Thompson 1994; Black- 
man & Yi 1998; Wheeler et al. 2000; Thompson, Chang, 
& Quataert 2004, hereafter TCQ04). Short-duration GRBs 
may result from black hole accretion-powered outflows fol- 
lowing the tidal disruption and coalescence of compact bina- 
ries (Paczyriski 1986, 1991; Eichler et al. 1989; Narayan et 
al. 1992; Ruffert et al. 1997; Janka et al. 1999). 

A unique property of both the surface of a newly-formed 
neutron star (a "proto-neutron star" or PNS; Burrows & 
Lattimer 1986) and in some cases the midplane of hyper- 
accreting disks is a significant excess of free neutrons (elec- 
tron fraction Y e <C 0.5), resulting from /^-equilibrium un- 

1 Astronomy Department and Theoretical Astrophysics Center, 601 
Campbell Hall, Berkeley, CA 94720; bmetzger@astro.berkeley.edu, 
eliot@astro.berkeley.edu 

- Department of Astrophysical Sciences, Peyton Hall-Ivy Lane, Princeton 
University, Princeton, NJ 08544; thomp@astro.princeton.edu 

3 Department of Physics, 366 LeConte Hall, University of California, 
Berkeley, CA 94720 

4 Lyman Spitzer Jr. Fellow 



der electron degeneracy (Pruet et al. 2003; Beloborodov 
2003a, hereafter B03a). 5 While degeneracy at a PNS's neu- 
trinosphere is assured, hyper-accreting disks are viscously 
heated and only possess a neutron excess in places where 
they are sufficiently dense and efficiently neutrino-cooled. If 
present at all, these neutrino-dominated accretion flows (or 
NDAFs) are thus geometrically-thin and confined to small 
radii in the disk (Popham et al. 1999). Recent neutrino-cooled 
a-disk calculations show that an NDAF only forms outside 
the last stable orbit for steady-state mass accretion above a 
critical "ignition" rate, given by 

M ig n ~ 0.07(0.02)A4 /3 ao / I 3 M Q s" 1 (1) 

for black hole mass M = 3Mt,Mq and spin a = 0(0.95), where 
a = O.lao.i is the disk viscosity parameter (Chen & Be- 
loborodov 2007, hereafter CB07). Disk mass accretion rates 
(Md) greater than M; gn , and thus neutron-rich NDAFs, are 
plausible in both collapsar and binary merger scenarios (e.g., 
Popham etal. 1999). 

For lower accretion rates or larger radii than charac- 
terize NDAFs, cooling is inefficient and the accretion is 
geometrically-thick and quasi-virial, forming an advection 
dominated accretion flow (Narayan & Yi 1994; Narayan et 
al. 2001). Although NDAFs must come into /3-equilibrium 
before accreting (B03a), this is not necessarily the case for 
thick disks; the neutron content of a thick disk may therefore 

5 In this paper we define a neutron excess as a neutron-to-proton ratio 
n/p = (\ -Y e )/Y e > 1, where we assume free nucleons and Y e is the proton- 
to-baryon ratio or "electron fraction." 
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depend on the matter that feeds it. We discuss the likely neu- 
tron content of thick disk winds more in £]5.5| and 8 16.31 

If the neutron excess that characterizes proto-magnetars and 
NDAFs can be preserved to large radii in the outflows that 
they power (where Y e obtains its asymptotic value Y"\ see 
Table Q] for definitions of commonly used variables), observ- 
able consequences may result. The dynamics of the GRB jet 
may be altered by the presence of a dominant neutron compo- 
nent (Derishev, Kocharovsky, and Kocharovsky 1999), which 
at proton-neutron decoupling could result in the emission of 
multi-GeV photons and neutrinos (Bahcall & Meszaros 2000; 
Meszaros & Rees 2000; Razzaque & Meszaros 2006a) and 
cause a significant reheating of the outflow (Rossi et al. 2006). 
Furthermore, a large neutron-to-proton ratio n/p may re- 
duce the fireball's asymptotic baryon contamination (Fuller et 
al. 2000; Vlahakis et al. 2003), contribute an additional com- 
ponent to the GRB light curve (Pruet & Dalai 2002; Bulik, 
Sikera, & Moderski 2002; Fan & Wei 2004; Fan et al. 2005; 
Dermer & Atoyan 2006) and a unique beta-decay signature 
(Razzaque & Meszaros 2006b), alter the subsequent afterglow 
emission (Beloborodov 2003b), and affect the fireball's nucle- 
osynthetic yield (Lemoine 2002; Pruet et al. 2002; B03a). Al- 
though not all of these consequences strictly require a neutron 
excess, most are far more conspicuous when n/p^> 1, partly 
because only excess neutrons will ultimately remain free if 
a-particle formation is efficient (e.g., Lemoine 2002). Iden- 
tification or strong upper limits on any of these effects would 
teach us much about the composition of GRB outflows. 

In this paper we examine the processes that shape the neu- 
tron content of outflows from GRB progenitors, motivated 
by the promise that the distinctive neutron-rich GRB signa- 
ture holds as a tool for constraining the elusive central en- 
gine, whose properties are masked by the outflow's otherwise 
rather generic dynamical evolution (e.g., the "fireball" model; 
Rees & Meszaros 1992). In particular, we focus on NDAFs 
and proto-magnetars rather than thick accretion disks, not be- 
cause these models are necessarily favored to produce GRBs, 
but rather because for neutron-rich central engines an asymp- 
totically neutron-rich outflow is plausible a priori. The goal 
of our analysis is to determine the conditions under which 
and degree to which these neutron-rich central engines can 
produce equally neutron-rich outflows. 

Determining the asymptotic neutron content of winds 
driven from the neutron-rich base of PNSs and NDAFs is non- 
trivial because the neutron fraction will evolve due to weak in- 
teractions under the comparatively nondegenerate conditions 
that characterize scales immediately larger than that of the 
central engine. In ^1 . ll we discuss the relevant processes that 
may "deneutronize" the outflow, driving n/p back to ~ 1. 
Indeed, in £|2] we show that thermally-driven outflows from 
PNSs and NDAFs are generally deneutronized by, if nothing 
else, electron neutrino absorption. Although thermally-driven 
winds of this kind possess little or no neutron excess, they 
also cannot produce GRBs because they do not reach ultra- 
relativistic speeds (Lorentz factor T ^S> 1). Winds driven di- 
rectly from the surface of PNSs or the midplane of NDAFs re- 
quire an energetically-dominant Poynting flux to reach large 
r because of the significant mass-loss driven by viscous and 
neutrino heating (e.g., Levinson & Eichler 1993); this re- 
quires rapid rotation and a strong magnetic field. If, through 
such enhanced magnetic acceleration matter is advected from 
the PNS surface or the NDAF midplane sufficiently rapidly 



to avoid deneutronization, the outflow will retain the large 
neutron-to-proton ratio that characterizes its degenerate base. 
The magnetocentrifugal acceleration required to maintain low 
Y e (high n/p), however, also enhances the wind's mass-loss 
rate over its purely thermally-driven value (TCQ04). This 
raises the question of whether simultaneously neutron-rich 
and ultra-relativistic outflows are possible under physically 
realizable conditions. 

To address these issues quantitatively we have calculated 
the neutron content of magnetically-driven winds from proto- 
magnetars and hyper-accreting NDAF disks by solving the 
equations of one-dimensional neutrino-heated magnetohydro- 
dynamics (MHD). We have studied in detail the effects of 
magnetic fields and rotation on PNS winds in a previous work, 
assuming a constant Y e (Metzger, Thompson, & Quataert 
2007, hereafter MTQ07). In <g]we include the evolution of 
Y e in these calculations in order to determine the asymptotic 
electron fraction Y" from proto-magnetar outflows. In <2]we 
adapt our calculations to the NDAF context by following out- 
flow from the accretion disk midplane for several flux tube an- 
gles, employing the a-disk NDAF models of CB07 as bound- 
ary conditions. We present a summary of our results in fj6j 
including a discussion of the prospect for neutron-rich out- 
flows from central engines of both long and short-duration 
GRBs. Our analysis and conclusions can be summarized in 
the broadest terms as follows: the conditions for n/p^$> 1 are 
simultaneously the conditions for short advection timescale, 
large mass-loss rate, and low asymptotic Lorentz factor. Thus 
only under very restrictive conditions do both n/p ^> 1 and 
high r obtain. The most promising possibilities appear to be 
outflows from sub-millisecond proto-magnetars and from ge- 
ometrically thick disks with Mo < M; gn (see Table|4|. 

1.1. Deneutronizing Processes 

Despite the presence of neutron-rich material at the base of 
NDAF and PNS outflows, this neutron excess may not be pre- 
served. Because conditions above the PNS surface and the 
accretion disk midplane are typically nondegenerate, equilib- 
rium between the pair-capture reactions 

e~ + p-*n + v e , (2) 

e + + n — > p + v e (3) 

favor Y e > 0.5 in the potentially pair-rich atmosphere through 
which the wind must accelerate. Furthermore, in the presence 
of an electron neutrino(antineutrino) energy density u ye {u^ e ) 
the inverse, neutrino absorption reactions 

v e +n— >e~+p, (4) 

i> e +p—>e + +n, (5) 

which dominate pair-captures abruptly above the outflow's 
launching surface, drive the system toward an asymptotic 
electron fraction given by 

Y »=( u gg. (^)-2A + 1.2AV( £Pf ) y' 
e ~ \ u Ve (e Ue )+2A + l.2Ay(e Ve )J ' W 

where A = 1 .293 MeV is the neutron-proton mass difference, 
(e„J((e p J) is the mean electron neutrino(antineutrino) energy 
(Qian et al. 1993; Qian & Woosley 1996, hereafter QW96), 
and the superscript v denotes that the electron fraction given 
by equation © is set solely by the properties of the local neu- 
trino radiation field. 
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Equation © shows that when the electron neutrino and 
antineutrino fluxes are comparable and have a similar spec- 
trum, as is generically the case for NDAFs and PNSs during 
the latter's Kelvin-Helmholtz cooling phase (Burrows & Lat- 
timer 1986), equilibrium between neutrino absorptions also 
favors a comparatively neutron-poor state (Y e ~ 0.5). 6 Such a 
deneutronizing luminosity of neutrinos must be present, self- 
consistently, for matter to cool to the dense, degenerate condi- 
tions required for low Y e in the first place. The total neutrino 
luminosity from a neutron-rich NDAF, for instance, must ex- 
ceed 

L„ )ign = r]M iga c 2 w 5 x lO^M^o^ 3 ergs s" 1 , (7) 

where r) w 0.04(0.15) for a = 0(0.95) (CB07) and we have 
used equation (HJ for Mj gn . Similarly, detailed numeri- 
cal calculations show that a cooling PNS's electron neu- 
trino/antineutrino luminosity is approximately given by 

L u P3 lO^Cf/lsr'ergss -1 (8) 

from a time t ss 1 s after core bounce until the end of the 
Kelvin-Helmholtz epoch at t = tkh ~ 10— 100 s (see, e.g., 
Pons et al. 1999, Fig. 14). 

To contrast the large neutron fraction in an NDAF's mid- 
plane with the much lower value favored in equilibrium with 
the NDAF's neutrino flux, Figure Q] shows the midplane elec- 
tron fraction Yf (dashed line) and the equilibrium electron 
fraction set by neutrino absorption (solid line) as a func- 
tion of disk cylindrical radius Ro (in units of gravitational radii 
R g = GM/c 2 ) for a steady-state NDAF solution taken from 
CB07 with a = 0.03, M = 3M Q , a = 0, and M D = 0.2M© s _1 . 
The local neutrino energy densities (u Ue ,u„ e ) and mean ener- 
gies ((ei/ e ),(ep e )) used to calculate from equation (O were 
obtained by integrating the total flux incident on a given po- 
sition just above the disk midplane at radius Ro, where Y" 
for an outflow launched near Rq is set. The disk neutrino 
emission is assumed to originate from axisymmetric annuli 
of negligible vertical height with a radial structure taken from 
CB07's one-dimensional, height-integrated calculations. We 
assume that all relevant lines of site are optically thin in cal- 
culating Y"\ this is a good approximation because the verti- 
cal neutrino optical depth through the disk is < 1 at all radii 
and because the disc scale height increases more rapidly than 
oc Rq so that the outer disk's atmosphere has an unobstructed 
view of the interior flow. Differential gravitational redshifts 
between emission and absorption radii are taken into account, 
although geodesic bending is ignored. 

Figure Q] shows that the disk midplane is very neutron-rich 
for R < 30- lOORg (the NDAF portion of the disk), reaching 
a neutron-to-proton ratio as large as ~ 30 at small radii. How- 
ever, FigureQ]also shows that Y" ~ 0.5 at all radii, so that if 
the outflow comes into equilibrium with the disk's neutrino 
luminosity it will be driven back to a relatively neutron-poor 
state (n/p ~ 1). Also note that F e " > 0.485 at all radii, al- 
lowing possible 56 Ni synthesis in a disk wind, again should 
the nucleons come into equilibrium with the neutrino flux 

6 During deleptonization, efficiently neutrino-cooled central engines re- 
lease slightly more u e 's than P e 's; in addition, the u e and v e spectra (and thus 
mean energies) also differ slightly due to the difference between the mean e~ 
and e + energies, the e~ and e + capture cross sections, and, in the neutrino op- 
tically thick case, the v e and v e neutrinosphere temperatures and geometries. 
Although Y" < 0.5 is possible in some cases, these relatively modest effects 
are unlikely toyield Y% < 0.5. 



(Pruet et al. 2004). 7 Since the outer disk is radiatively in- 
efficient and therefore particularly prone to large mass out- 
flows (Blandford & Begelman 1999), disk wind-aided stel- 
lar explosions provide one way to produce optically luminous 
SNe in collapsar models for LGRBs (MacFadyen & Woosley 
1999), where Ni masses up to ~ 0.5Mq have been inferred 
(e.g., GRB980425/SN1998bw; Iwamoto et al. 1998; Woosley, 
Eastman, & Schmidt 1999). 

2. THERMALLY-DRIVEN WINDS 

In spite of nondegenerate pair captures and neutrino absorp- 
tions (eqs. |2]-||5)), the high neutron fraction in NDAFs and 
proto-magnetar outflows will be preserved if deneutronization 
proceeds slower than a fluid element's advection from the sur- 
face. In evaluating this possibility in the case of NDAFs, 
previous investigations have assumed that the disk is turbu- 
lent and that the relevant outflow rate is of the order of the 
turnover frequency of a typical turbulent eddy (~ f^, the 
Keplerian rotation rate), concluding that low Y e is preserved 
under most conditions appropriate to hyper- accreting disks 
(Pruet et al. 2003; B03a). Specifically, Pruet et al. (2003) 
argue that because NDAFs are generally dense, degenerate, 
and gas pressure-dominated, if a fluid element is carried out 
of the midplane adiabatically in an eddy turnover time, pair 
creation is somewhat suppressed and low Y e is preserved (i.e., 
Y" sa F C D obtains). B03a also compares JIk to the neutrino 
absorption rates (eqs. J!)-[|5)), reaching a similar conclusion 
to Pruet et al. (2003) except for very high disk neutrino lu- 
minosities (L v > 10 53 erg s" 1 ), corresponding to Mq > 1M 
s -1 . This leaves a wide range of astrophysically plausible ac- 
cretion rates (M; gn < Mp < 1M Q s" 1 ) for which outflow, if it 
escapes the disk in an eddy turnover time, would be asymp- 
totically neutron-rich. 

Several arguments can, however, be raised against the con- 
clusion that low Y e is preserved in NDAF outflows. Although 
some degree of turbulence and turbulent mixing must ac- 
company the accretion torque (e.g., via the MRI; Balbus & 
Hawley 1998), its scale and efficiency are unclear. 8 Even if 
present, large-scale turbulence is not likely to drive outflow 
from a thin disk; thus, an eddy turnover time is not the rele- 
vant timescale to compare to the weak interaction rates. Out- 
flow from the surface of a thin accretion disk must be heated, 
passing through a sonic point as it accelerates out of the black 
hole's potential. Thus, the advection time of a self-consistent, 
viscous or neutrino-heated wind sets the residence time for a 
fluid element in regions of potential deneutronization. This 
timescale is generally much longer than . 

Indeed, the neutrino-driven PNS winds that accompany the 
Kelvin-Helmholtz cooling of slowly rotating, non-magnetized 
PNSs generically come into equilibrium with the neutrino 

7 Although we find that Y" > 0.485 at all disk radii using mean neutrino 
energies from CB07's height-integrated disk calculations, precisely whether 
Y" > 0.485 or Y" < 0.485 is difficult to determine with confidence because 
the neutrino spectra are sensitive to the disk's vertical temperature profile, 
which is theoretically uncertain, and to neutrino transport if the disk midplane 
is neutrino optically thick (Sawyer 2003). If viscous heating is important in 
the wind, additional entropy deposition may drive Y" > 0.485, allowing 56 Ni 
production even if Yj? < 0.485. 

8 NDAFs cool efficiently and are not unstable to radial convection, except 
perhaps at small radii in the case of neutrino optically thick accretion (Lee 
et al. 2005); it is also not clear that NDAFs are vertically-convective (Hi- 
rose, Krolik, & Stone 2006), as this depends on the unknown vertical viscous 
dissipation profile (see the discussion in Blaes 2007). 
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FIG. 1. — [LEFT AXIS] Electron fraction in neutrino absorption equilibrium Y" (eq. (§); solid line) as a function of the wind launching radius So (in units 
of gravitational radii R g = GM/c 2 ) calculated from the thin, height-integrated a-disk accretion model of Chen & Beloborodov (2007) with a = 0.03, black hole 
mass M = 3Mq , black hole spin a = 0, and accretion rate Md = 0.2 Mq s -1 . Also shown is the midplane electron fraction Y® (dashed line) taken directly from 
CB07's calculation. Note that the disk midplane is very neutron-rich inside Rq ~ 30R g , with a neutron-to-proton ratio n/p = (\—Y e )/Y e exceeding 10. However, 
if outflow driven from the disk comes into equilibrium with the disk's neutrino radiation then the asymptotic neutron-to-proton ratio is driven back to ~ 1 because 
Y" ~ 0.5 at all radii. Also note that neutrino equilibrium favors Y" > 0.485 at large radii, a requisite condition for producing 56 Ni in neutrino-driven outflows. 
[RIGHT AXIS] The dot-dash(triple-dot-dash) line shows the mean electron neutrino(antineutrino) energy incident on the base of outflows driven from radius 
Rq. Also shown (dotted line) is the binding energy Eg of a nucleon in a Keplerian thin disk, where the black hole's gravitational potential is assumed to be 
Newtonian for simplicity. The disk is the most neutron-rich at radii where Eb > (tv e ), {^u e ) and thus where more than one neutrino absorption is required to 
unbind a nucleon. This shows that purely neutrino-driven winds from small radii in the disk must come into equilibrium with neutrino absorptions and will thus 
obtain Y" ~ Y" ~ 0.5 (n/p ~ 1). At larger radii (Rq ~ 50 — 200R g ), Eb < {tv e ), and so deneutralization of a wind by neutrino absorption is less likely. 



absorption rates, obtaining Y" = Y" (QW96; Thompson et 
al. 2001, hereafter T01). To see that the neutron fraction of a 
purely neutrino-driven NDAF wind must also come into equi- 
librium with the neutrino rates we adapt an argument first pre- 
sented by QW96 in the PNS context. First, note that a typical 
nucleon with mass m n launched in a wind from radius Rq of 
a thin Keplerian disk requires an energy Eb — GMm n /2Ro « 
80(Ro/6R g )~ 1 MeV to escape the black hole's gravitational 
potential, where we have assumed Newtonian gravity for sim- 
plicity. In addition to Yf and F e ", in Figure Q] we show for 
comparison Eg and the mean electron neutrino and antineu- 
trino energies used to calculate Y" . The mean neutrino ener- 
gies vary only weakly with radius because the neutrino flux 
released from an NDAF's inner disk is relatively hot (Popham 
et al. 1999; CB07) and dominates the neutrino heating at all 
outflow radii (relative to the neutrino emission from larger 
radii in the disk). By balancing viscous heating and opti- 
cally thin neutrino cooling at the radius R p ~ 1 - 10R g where 
most of the disk's emission originates, we find that an NDAF's 
mean neutrino energy incident on any disk radius is approxi- 



mately given by 



'5.17b 



1/5, 



-1/5, 



(R p /6R g T 3/10 MeV, (9) 



where To is the disk's midplane temperature at R p . Equa- 
tion (0, which agrees well with our calculation of {e u ) in 
Figure [1] is valid so long as the inner disk is optically thin 
to neutrinos (t„ < 1), which remains true for accretion rates 
below 0.7(0.05)ao iM Q s" 1 for electron neutrinos and w 
2.0(0. 14)a .i^ Q s" 1 for electron antineutrinos for a = 0(0.95) 
and M = 3M Q (CB07). For larger M D the disk is opaque to 
neutrinos (T Ue ,T^ e > 1) and the mean neutrino energy is ap- 
proximately given by its Fermi-Dirac blackbody value 

(e„) w 3.27; » 16M3 1/2 M^o.i(^ P /6^r 1/2 MeV, (10) 

where Md = 0. IMd.o.iMq s_1 an d T v is the temperature at the 
disk's neutrinosphere. 

Despite the relatively high mean neutrino energies shown in 
Figure[T]and implied by equations d9b-([T0b. (e„) is still much 
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less than Eb(^ R p ) because NDAFs are efficiently neutrino- 
cooled. Thus, each nucleon in a thermally-driven outflow 
originating from an inner NDAF must absorb several neutri- 
nos in the process of being unbound from the black hole's 
potential. Through these multiple absorptions the wind's elec- 
tron fraction is unavoidably brought into equilibrium with the 
neutrino flux at Y" (QW96), a consequence of the fact that 
neutrino absorptions (eqs. H-lO) both dominate the wind 
heating and determine Y" ~ Y" , 

This conclusion may hold even for winds driven from near 
the NDAF's relatively loosely -bound outer edge. NDAFs only 
exist interior to an "ignition" radius R lga , which we estimate 
as the location at which optically thin cooling by nondegener- 
ate pair capture on free nuclei (eqs. EJ-JS) balances viscous 
heating for a thick disk: 9 



Ru 



(11) 



For the disk parameters associated with the solution in Fig- 
ure [J (q = 0.03, M = 3M Q ,M D = 0.2M Q s" 1 > M ign ), equa- 
tion (fTTT i gives Ri gn w 100/? g , corresponding to the radius in- 
terior to which the disk possesses a significant neutron excess 
(Yf <C 0.5). By requiring that Ri gn exceed the radius of the 
innermost stable circular orbit, /?; sco ~ 6(l)R g for black hole 
spin a = 0(0.95), one recovers the numerically-determined 
value for M{„ n given in equation (HJ to reasonable accuracy. 

Using equations (0 and ([Tol l, the ratio of the binding en- 
ergy at R lgn to the mean neutrino energy from the inner disk is 
approximately given by 



-^D.O.l 3 "0.1 I 



3/10 



< 1 



o.i 



6R< 



1/2 



■r v >l 



(12) 

Equation dT2b shows that EsiRigx) > (e„) for a > 0.1 and 
for most physical values of Md; thus, winds thermally-driven 
from even Rq ~ R ign would have Y" w Y" ~ 0.5. Although 
equation (fT2l implies that winds driven from the outer por- 
tions of high-Mo, low-a NDAFs may remain neutron-rich, the 
astrophysical situations most likely to result in extended disks, 
collapsars and He core-black hole mergers (Fryer & Woosley 
1998), are also likely to have lower values of M D . Further- 
more, NDAFs become unstable to self-gravity (Toomre Q < 1) 

atR lga for Mo > 2a^( 7 Ml^ 1 M & s" 1 , comparable to the accre- 
tion rates for which EB(Ri gn ) < (e^)- The dynamics at these 
radii is not well understood so it is difficult to draw defini- 
tive conclusions concerning the neutron content of outflows 
launched from large radii in high-Mo, low-a disks. 

In addition to neutrino heating, viscous heating may ther- 
mally drive outflows from the surface of NDAFs. Pruet et 
al. (2004) have shown, in the context of a wind driven by vis- 
cous heating supplied through a simple a prescription, that 
entropy added to the outflow lifts electron degeneracy be- 
fore the outflow falls out of /3-equilibrium. This drives n/p 

9 Equation jilt overestimates R Im for large R[„ n (i.e., large Mo and low a) 
because cooling via electron capture requires a threshold energy A— nteC 2 ~ 1 
MeV that exceeds the disk temperature at large radii. Furthermore, NSE 
favors composite nuclei for Rq > R a ~ 100R g ; nuclear disintegration is a 
significant heat source around So ~ Ra and we have overestimated cooling 
for Rq > R a because pair capture rates are lower on composite nuclei. 



to a value near unity via nondegenerate pair captures, even 
without the aid of neutrino absorption. However, the abil- 
ity of viscous heating to generically drive outflows from thin 
disks is unknown, as this depends in detail on where energy 
is deposited, a major uncertainty in current thin disk theory 
(Blaes 2007). Recent radiation MHD thin disk simulations 
suggest that very little energy dissipation occurs in the disk 
corona (Turner 2004; Hirose, Krolick, & Stone 2006; Blaes 
et al. 2006b; Rrolik et al. 2007), where a wind would most 
likely be launched. Regardless of potential viscous entropy 
contributions, the existence of the neutrino flux discussed in 
the preceding paragraphs is assured. We thus conclude that 
purely thermally-driven winds from the vicinity of compact 
objects, even those with a significant surface neutron excess, 
are unlikely to be neutron-rich asymptotically. 

3. MAGNETICALLY-DRIVEN WINDS 

Winds driven directly from the surfaces of NDAFs and 
PNSs solely by viscous or neutrino heating are inherently 
mass-loaded and thus make poor candidates for producing 
GRBs. The high energy-to-baryon ratio required of ultra- 
relativistic GRB outflows (asymptotic wind Lorentz factor 
T > 10- 100; e.g., Lithwick & Sari 2001) can, however, be 
achieved if rapid rotation and a strong poloidal magnetic field 
supply an energetically dominant Poynting flux. The ratio 
of Poynting flux to kinetic energy flux at the light cylinder 
/?l = c/f2 is given by (for a > 1 ) 



B 2 



(13) 



where B is the magnetic field strength, p is the rest mass den- 
sity, SI is the rotation rate of the central star or disk, and 
we have assumed that the outflow is moving only mildly 
relativistically at R^. If the magnetic energy is fully con- 
verted into the kinetic energy of bulk motion, either directly 
or through thermalization and subsequent thermal or magnetic 
pressure-driven expansion (e.g., Spruit & Drenkhahn 2002), 
then r ~ cr; Poynting-flux dominated GRB outflows therefore 
require cr > 10- 100. 

In the case of magnetically-driven, high-cr NDAF or PNS 
winds, matter may be advected from the wind's base suffi- 
ciently quickly to remain effectively adiabatic, with its initial 
degeneracy never lifted due to insufficient heating. The pair- 
capture reactions (eqs. EJ-O) then continue to favor low Y e 
above the PNS surface or disk midplane. In the PNS case this 
ineffective heating becomes manifest as an exponential drop 
in the asymptotic wind entropy with increasing il, which oc- 
curs for rotation periods P = 2ir/fl < 2-3 ms (see eq. [37] 
of MTQ07) in the presence of a sufficiently strong surface 
dipole magnetic field (Z?J P ~ 10 14 -10 15 G); note that the large 
fields and rapid rotation required are similar to those required 
for proto-magnetars to produce LGRBs in the first place. It 
is therefore likely that GRB jets from proto-magnetars or 
NDAFs are not completely deneutronized by pair captures. 

Even with degeneracy intact and pair production sup- 
pressed, to remain neutron-rich a GRB-producing wind must 
advect material from its base sufficiently rapidly to over- 
come neutrino absorptions. This requires significant mag- 
netocentrifugal support in the outflow's inner, hydrostatic at- 
mosphere, where Y" is set. For the reasons discussed in SjH 
magnetocentrifugal forces acting on a wind launched from 
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a radius R must contribute a factor ~ GMm n / R(e v ) ^> 1 
more energy than neutrino heating in unbinding the outflow 
to avoid deneutronization. Inevitably, such acceleration in 
the wind's subsonic region leads to significantly enhanced 
mass-loss (TCQ04; MTQ07). Thus, the very conditions 
required to preserve low Y e in a proto — magnetar or NDAF 
wind threaten to simultaneously over — pollute the outflow 
with baryons, reducing a and stif lying the outflow's GRB 
potential. 

In the following sections we address some of these issues by 
calculating the neutron content of magnetized proto-magnetar 
and NDAF winds. In particular, we consider the conditions 
under which outflows from proto-magnetars and NDAFs 
(© can remain neutron-rich while simultaneously maintain- 
ing sufficiently low mass-loading to remain plausible GRB 
central engines. 

4. PROTO-MAGNETAR WINDS 

4. 1 . Evolution Equations and Numerical Procedure 

As in our previous work (MTQ07), we calculate the struc- 
ture of rapidly rotating PNS winds by solving the equations 
of one-dimensional neutrino-heated, ideal MHD in the equa- 
torial plane of the PNS. Using the time-dependent "inhomo- 
geneous" 2N-RK3 scheme described in Brandenburg (2003), 
we solve for the wind density p, temperature T, radial veloc- 
ity v r , azimuthal magnetic field Bq, and azimuthal velocity 
as a function of radius r according to equations (2)-(5) and 
equation (9) in MTQ07. In this work we extend our previ- 
ous calculations by simultaneously solving for the electron 
fraction Y e , thus determining the wind's asymptotic neutron 
abundance. 

The electron fraction Y e evolves as the wind emerges off the 
PNS surface due to weak interactions according to: 

d 

+ \<rp^n»,), (14) 

where d/dt = d/dt + v r (d/dr) and the A's are the weak in- 
teraction rates; we take the pair capture rates (eqs. HJ-O) 
and the neutrino capture rates (eqs. ID-El) from B03a (ne- 
glecting proton/neutron recoil to good approximation; Stru- 
mia & Vissani 2003). We have included the full effects 
of electron degeneracy, including the outgoing electron and 
positron blocking factors in calculating the neutrino capture 
rates. In calculating all of the weak interaction rates we as- 
sume that the electrons and positrons are relativistic. This is 
a good approximation for neutrino capture because the av- 
erage neutrino(antineutrino) energy (and hence the resultant 
electron(positron) kinetic energy) always far exceeds the elec- 
tron rest mass. We set the pair capture rates equal to zero 
for T < 0.5 MeV in order to artificially account for the dis- 
appearance of pairs; the evolution of Y e is not sensitive to 
this cutoff because for T ~ 0.5 MeV pair capture is always 
dominated by neutrino absorption. We neglect the effects that 
a-particle formation has on the evolution of Y e ; this is a rea- 
sonable approximation because for cases in which the formed 
a— particle fraction is the most significant (i.e., Y" ~ 0.5) and 
is thus likely to have its greatest effect on the evolution of Y e , 
Y" almost obtains by the radius where most a-particles form. 
We also neglect the effects that magnetic fields have on the 
electron and positron distribution functions and, hence, on the 



interaction rates and equation of state, although these effects 
become important for the largest field strengths that we con- 
sider (B = 10 16 G) and should be included in more detailed 
work (e.g., Lai & Qian 1998; Duan & Qian 2004). Lastly, 
we neglect the small effect that general relativity (GR) has on 
the evolution of Y e ; in non-rotating PNS winds GR slightly 
increases Y" (Fuller & Qian 1996), primarily due to neutrino 
gravitational redshifts and the deeper gravitational potential. 

Because we are interested in the wind structure along open 
magnetic flux, we assume a monopole radial field structure, 
B r = B v (R v /r) 2 , where B v is the surface magnetic field, and 
R v is the neutrinosphere radius. Since Y" is set near the PNS 
surface, Y e is more sensitive to the surface field strength than 
to the field's precise radial scaling, and so our results are likely 
relatively insensitive to our monopole assumption. 

For asymptotically relativistic PNS outflows, using the con- 
served magnetic flux <£>b = B V R 2 = B(Ri)R^ and a spherically 
symmetric mass flux M = 4-npv r r 2 = AircR^piRi} we have that 
the magnetization from equation ( fT3l is given by 

a = B 2 Xn 2 /Mc 3 (15) 

Although for asymptotically non-relativistic outflows (a < 
1), a as defined by equation <TT3T > is no longer the ratio of 
Poynting-to-kinetic energy flux at R^, a remains a useful di- 
mensionless quantity for characterizing the importance of the 
magnetic field in accelerating the outflow. 

With the exception of the evolution of Y e and its effects on 
the other thermodynamic wind variables, our boundary con- 
ditions and microphysics are similar to those used in MTQ07; 
thus, we review only the most essential aspects of these 
here. Our neutrino heating(cooling) rates include charged- 
current neutrino absorption(pair capture) (eqs. J2]-||5)), neu- 
trino(pair) annihilation, and inelastic neutrino-lepton scatter- 
ing (QW96; T01), corrected for solid angle and redshift ef- 
fects (Salmonson & Wilson 1999; T01). We index stages 
of the PNS thermal evolution in terms of the electron an- 
tineutrino luminosity L Pe . We scale all other neutrino lu- 
minosities (L„ f ,L„ ,L P ,L Vt , and L Pt ) as in TCQ04: L„ e = 
LvJ 4.3 = 1.08L„ , where p denotes each of the other four 
neutrino/antineutrino species. Note that the total neutrino lu- 
minosity is then L v ~ 4.6L Pf . Following T01, all first energy 
moments at the neutrinosphere ((e„) = (E 2 ) / (E v ), where E v 
is the neutrino energy) were scaled with luminosity as (e„) oc 

L X J\ anchoring {(e Ve ), (e Pe ), <e„ M }} at {11, 14,23 }MeV for 
Lv e .5\ = 8, where Lp f ,5i is the electron antineutrino luminos- 
ity in units of 10 51 ergs s . Higher energy moments neces- 
sary for the heating calculations ((e" ), (e^ ), etc.) are related 
to the first through appropriate integrals over the assumed 
Fermi-Dirac surface distribution. With these scalings we note 
that electron antineutrino luminosities Lp„5i = {8,3.5, 1} re- 
sult in a neutrino-driven asymptotic electron fraction = 
{0.48,0.50,0.53} according to equation ©; neutrino absorp- 
tions therefore favor n/p ~ 1 throughout the PNS's Kelvin- 
Helmholtz cooling epoch (T01). Although we have assumed 
neutrino spectra based on current calculations of non-rotating 
PNSs, for rotating PNSs the mean neutrino energies along 
the equator may be lower due to gravity-darkening, which 
would increase Y" from the non-rotating value (e.g., Fryer & 
Heger 2000; Thompson, Quataert, & Burrows 2005; Dessart 
et al. 2006). 

We assume Newtonian gravity for a fixed central PNS mass 
M = 1 AMq and a neutrinosphere radius R„ = 10 km which is 



characteristic of the PNS's final, cooled state. Although /?„ is 
probably larger than 10 km for the first few seconds following 
the launch of the SN shock (e.g., Buras et al. 2003), neutron- 
rich GRB-producing outflows, which are the focus of this pa- 
per, are only possible at relatively late times, once the PNS has 
fully contracted to its most rapidly rotating state (see Ej4.4| >. 
Our code is non-relativistic but we still calculate flows with 
a > 1, which accelerate to relativistic speeds outside the light 
cylinder, because M and Y e are set very close to the PNS sur- 
face, where the wind still moves non-relativistically. We set 
the neutrino sphere density of the wind p v so that the neutrino 
optical depth to infinity is t v — \\ Pv ranges from ~ 10 12 g 
cm" 3 for high luminosity, rapidly rotating solutions to > 10 13 
g cm -3 for our lowest luminosity solutions (see Table [2]). In 
general, we find that Y° is relatively insensitive to p v . The 
electron fraction at the neutrinosphere F c ° is chosen to ensure 
that equilibrium between the weak interaction rates in equa- 
tion ( TBI is established at R v . 

The value of the magnetic flux <£>b = B V R\, stellar rotation 
rate f2, and neutrino luminosity L v uniquely identify a wind 
solution. By letting the wind come into steady-state, we ob- 
tain eigenvalues for the mass-loss rate M (normalized to 4ir 
sr), specific angular momentum loss rate C, and Bernoulli 
integral B (see eqs. [6], [7], and [11] in MTQ07). Because 
our code is non-conservative, the radial conservation of these 
quantities is used to verify the code's accuracy. The total 
asymptotic energy lost in the wind is given by E a = B a M, 
where B a is the Bernoulli integral evaluated at large radii. 
The rate of angular momentum loss in the wind is given by 
Jw = CM = QR\M, where R& is the Alfven radius defined by 
B r {R A )/^4np(R A ) = v r {R A ). 

4.2. Numerical Results 

We have calculated the MHD structure of PNS winds for 
several combinations of surface monopole magnetic field 
strength B v , rotation rate ft, and electron antineutrino lumi- 
nosity Lp e = L Ve .51 x 10 51 ergs s" 1 in order to study the neutron 
fraction in proto-magnetar outflows. Some of these results are 
summarized in Table [2] 

Figure [2] shows our calculation of the asymptotic electron 
fraction Y" of proto-magnetar winds as a function of il at neu- 
trino luminosity L Ve >51 = 8 and 3.5 forB,, = 10 14 , 10 15 , and 10 16 
G. Figure [2] and Table [2] show that for slow rotation and low 
magnetic field strengths Y" obtains the neutrino absorption 
equilibrium value Y e v ~ 0.5 (eq. J6)), as is expected from stud- 
ies of slowly-rotating, non-magnetized PNSs (QW96; T01) 
and from the arguments given in S|2] As the rotation rate 
and magnetic field strength increase, Y" decreases because 
the material is advected sufficiently rapidly from the surface 
due to magnetocentrifugal slinging that Y e "freezes out" be- 
fore coming into equilibrium with neutrino absorptions. 

Figure |2] also shows that for a sufficiently strong magnetic 
field, B v > 10 14 - 10 15 G, Y e a (fl) no longer increases with in- 
creasing B v , saturating to a profile F e fl sat (fi) that we find is 
reasonably well-fit by a single empirical formula for all of the 
rotation rates and neutrino luminosities that we have consid- 
ered (Q, < 9000 s" 1 ; 1 < Lp e) 5i < 10, corresponding to times 
~ 1 - 10 s following bounce): 
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FIG. 2. — Asymptotic electron fraction Y" for magnetocentrifugally-driven 
PNS winds as a function of rotation rate Q at La e = 8 X 10 51 ergs s -1 and 
Lp e = 3.5 X 10 51 ergs s" 1 for B v = 1() 14 G, B v = 10 15 G, and B v = 10 16 
G. For slow rotation and weak magnetic fields, Y° approaches its neutrino 
absorption equilibrium value Y" ~ 0.5 (eq. (6)). For more rapid rotation 
and stronger magnetic fields, Y° is reduced because matter is advected from 
the PNS surface sufficiently rapidly that it falls out of /3-equilibrium before 
neutrino absorptions dominate degenerate pair captures (see Fig. f5J. For 
sufficiently large B v , Y" no longer decreases with increasing B v because the 
wind corotates past a few scale heights above the PNS surface, where Y" 
is set. A n ap proximate fit to our numerical results in this limit is given by 
equation 1161 . 

where fi n w 7800 s~\ Afi n « 2000 s" 1 , and Y e ° is the elec- 
tron fraction at the neutrinosphere (typically, Y® w 0.01-0.05; 
see Table |2). Equation ( fTo*l ) and Figure |2] show that proto- 
magnetars must have submillisecond (P < P n = 2n/D, n rs 0.8 
ms), near break-up, rotation in order to produce asymptoti- 
cally neutron-rich outflows; furthermore, because ASl n /f2 n <C 
1 the transition from n/p ^> 1 to n/p ~ 1 occurs over a very 
limited range in Q. Therefore, because a PNS possesses a 
rotational energy 



^rot — 



4x10 



52 



ergs, 



(17) 



(16) 



where / w (2/5)MR 2 u is the PNS moment of inertia, we con- 
clude that neutron-rich outflows from magnetar birth would 
require a GRB plus SN energy totaling >4 x 10 52 ergs. 

As an example of a neutron-rich wind solution, the left 
panel of Figure [3] shows the electron fraction Y e as a func- 
tion of radius just above the PNS surface for a solution with 
L Pe 51 = 8,B V = 10 16 G, and fi = 9000 s" 1 (rotation period 
P as 0.7 ms). This relatively high luminosity solution corre- 
sponds to a time wis after core bounce (see eq. (HI). Figure 
[3] also shows Y" and the equilibrium electron fraction T c eq ob- 
tained by setting the time derivative on the left hand side of 
equation ( TT4b equal to zero. The right panel of Figure[3]shows 
the individual rates in the wind that determine the evolution 
of Y e : electron absorption (eq. J2, dotted line), positron ab- 
sorption (eq. Q, dashed line), electron neutrino absorption 
(eq. HI, triple -dot -dashed line), electron antineutrino ab- 
sorption (eq. 0; solid line), and the wind's advection rate 
A a dv = v r J r (dot-dashed line). The slight difference between 
the neutrino and antineutrino absorption rates at large radii, 
although difficult to discern on the scale of this plot, implies 
that Y? ~ 0.48 for this solution. 

Figure [3] shows that the electron fraction at the base of the 
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FIG. 3. — Left: Electron fraction Y e (solid line) in the inner 6 km of a magnetically-driven PNS wind with M = 1.4Mq, R v = 10 km, B v = 10 16 G, and 
f2 = 9000 s -1 ; also shown is the local equilibrium electron fraction (long dashed line), obtained by setting the left hand side of equation )14t equal to zero. 
The electron fraction is very small at the PNS surface (Y e 0.05). As matter accelerates away from the surface, Y e m Yg as long as matter moves sufficiently 
slowly to maintain local equilibrium with the weak interaction rates. Because Y e is out of equilibrium by the time neutrino absorptions dominate pair captures, 
the asymptotic electron fraction Y" is <C 0.5, despite the fact that Y? q asymptotes to Y" ~ 0.48 at large radii. Right: Rates in the wind that determine the radial 
evolution of Y e . The dotted(short dashed) lines show the electron(positron) absorption rates (eqs. |5]-(3]). The electron neutrino and antineutrino absorption rates 
(eqs. (4)— (5)) are shown with a triple-dot-dash and solid line, respectively; the wind advection rate A a[ i v = v, /r is shown with a dot-dashed line. The asymptotic 
electron fraction obtained by the wind Y" ss 0.19 is appreciably lower than the neutrino absorption equilibrium value (Y" ~ 0.48) that obtains for slower rotating 
or non-magnetized winds. Lower Y° obtains because magnetocentrifugal acceleration is sufficiently strong that matter "freezes out" of /3-equilibrium at small 
radii (corresponding to A a( j v rising above the weak interaction rates at r ss 12 km), before the neutrino absorption rates begin to dominate the electron capture 
rate. Equilibrium still favors neutron-rich matter at freeze-out because matter is advected from the PNS surface sufficiently quickly that neutrino heating is 
suppressed and degeneracy is never completely lifted. 



outflow is very neutron-rich (Y e ° ps 0.05). This results from 
the equilibrium established between the electron and electron 
neutrino capture rates (dotted and triple- dot -dashed lines, 
respectively) under the degenerate surface conditions (p„ sw 
10 12 g cm" 3 , T„ ps 4.2 MeV < 14 MeV ps fi e , where (i e is 
the electron chemical potential at the neutrinosphere), which 
strongly suppresses the thermal positron population relative to 
electrons. Thus our value of F e ° is significantly larger than the 
pure pair-capture equilibrium value given in Figure 1 of B03a 
because we have included neutrino absorptions; for vanishing 
neutrino luminosity, however, the results of B03a Figure 1 
would obtain. Indeed, for solutions with lower L Ue = \3L Ve 
electron neutrino absorption becomes less important and Y® 
decreases (see Tabled. 

As matter accelerates off of the PNS surface, Figure [3] 
shows that Y e rises above Because the wind initially 
moves slowly relative to the relevant weak interaction rates, at 
small radii Y e remains in approximate equilibrium at the value 
Y e ps F(f q < Yg. However, as magnetocentrifugal slinging ac- 
celerates matter away from the PNS surface, A a dv = v r /r also 
rises rapidly, eclipsing the weak interaction rates by r ps 12 
km. At this point Y e "freezes out" at a value Y" ~ 0.19; this 
occurs before neutrino absorptions dominate the electron cap- 
ture rates (which is why Y e ps T t f q is still relatively low). 

4.3. Conditions for Neutron-Rich Outflows from 
Proto-Magnetars 

The mass-loss rate for purely thermal, neutrino-driven PNS 
winds M(a can be derived analytically by requiring that the en- 
ergy used to unbind a typical nucleon Eb = GMm n /R u be sup- 
plied entirely by the neutrino heating Q v = J q u dr/v r that 
a nucleon experiences in accelerating from the PNS surface 
to large radii (QW96), where q v is the net neutrino heating 
rate per baryon, which is dominated by neutrino absorption 



(eqs. flU-G)). By repeating this derivation but instead requir- 
ing that Q u < (e v ) (see Sj2]l we obtain the minimum mass-loss 
rate that must accompany a neutron excess over the neutrino- 
driven value 

M ^ GMm n , 

-T- > <f> n = ~ 20M L4 e Pf ,io " ^To, (18) 

M th Ru\Ev) 

where (QW96) 

Mth w lO- 6 4f 51 (e Pf , 1 o) 1O/ X.4^io 3 M s- 1 , (19) 

R u = 10R w km,M= lAM h4 M Q , and (e u ) sa (e p .) = 10(e p „io) 
MeV. In our numerical calculations we have assumed that 
L v oc (e„) 4 so that M t h L 5 J 2 . 

In Figure [4] we show M for the same wind solutions for 
which Y" is shown in Figure|2] Figure [4] shows that M is sub- 
stantially enhanced over its thermally-driven value M t h in the 
presence of a strong magnetic field and rapid rotation. This 
enhanced mass-loss occurs because a strong magnetic field 
forces the outflow to corotate above the PNS surface and rapid 
corotation brings the wind's sonic radius R s much closer to the 
surface, into the corotating region; this increases the hydro- 
static scale height and the wind density at /? s , thereby increas- 
ing M. Because the magnetic field, as opposed to neutrinos, is 
primary responsible for unbinding matter from a rapidly ro- 
tating proto-magnetar, the outflow can have M > M t h<pn and 
can thus remain neutron-rich. Indeed, Figure [4] shows that M 
increases exponentially with £1 and a comparison with Fig- 
ure |2] shows that Y" first noticeably decreases below once 
equation ( fl~8T > is satisfied. 

Figure|4]also shows that for modest magnetic field strengths 
M increases with increasing B„, but that for sufficiently large 
B v , M saturates to a value M c f in the "centrifugal limit" (e.g., 
Lamers & Cassinelli 1999). For il ps fl n (or, equivalently, 
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FIG. 4. — Mass loss rate M for the same solutions for which Y° is shown 
in Figure [2] Due to centrifugal slinging, M increases with B v and Q and 
is substantially larger than its thermally-driven, non-rotating value, which 
is M m K. 2 x 10r*(2 X 10 _5 )M Q s" 1 for = 8(3.5). By comparing M 

with Y" in Figure|2]note that winds with large M have low Y"\ in particular, 
M > (f> n M m ~ 15-20M C „ (eq. 1181 ) is required for Y" to noticeably decrease 
below its neutrino-driven value Y" . 

Y? w 0.25) we find that 

M cf (fU « 6 x 10- 2 (L Pf ,5i/8) 2 2 M Q s _1 , (20) 

which is over an order of magnitude greater than the mass-loss 
required, M th n «3x 10- 3 (L Pf , 5 i/8) 225 M o s~ l (eq. ED), for 
just a mild neutron excess relative to the neutrino-driven value 
(i.e., Y e a < Y»). 

Figures |2] and |U show that at fixed fi, M — > M cf and Y" — > 
ya.sat (16]) for similar magnetic field strengths. To un- 
derstand why this occurs note that the requirement for fully 
centrifugally-enhanced mass-loss is that the magnetic field 
must be sufficiently strong to enforce corotation beyond the 
radius where M is set; because a wind approximately coro- 
tates to its Alfven radius Ra and its mass-loss rate is set at its 
sonic radius R s the condition for M — > M c f is that Ra >R s - In 
analogy, because Y" obtains near the PNS surface (see Fig. [3), 
the condition for Y" — > Y"' sdt is that R A >R„. In the centrifu- 
gal limit the sonic radius for an equatorial wind is given by 
R s = (GM/fi 2 ) 1 / 3 (e.g., Lamers & Cassinelli 1999), so that 
R s (£l n ) « 13 km ~ R v . Thus the conditions for fully enhanced 
mass-loss and fully enhanced neutron-richness are approxi- 
mately the same because the sonic radius is close to the PNS 
surface for the near break-up rotation rates required to pro- 
duce neutron-rich outflow. 

Finally, although we have focused on equatorial proto- 
magnetar winds, the properties of the outflow will vary with 
latitude. In particular, M and Y" will decrease and increase 
with increasing angle 9 from the equator, respectively, 
approaching M± and Y" for nearly polar outflow (6 — > 90°) 
because of the reduced centrifugal acceleration along nearly 
vertical field lines (M and M t h here are per solid angle). 
Thus, it may seem plausible that outflows driven from 
intermediate latitudes could maintain moderately low Y e yet 
become significantly more relativistic than matter driven 
from near the equator. Even high-latitude outflows, however, 
require M » M± (eq. ED) to attain Y" < F e ", making such a 
scenario unlikely, as we will show explicitly through NDAF 



disk wind calculations in S|5] Furthermore, the limited solid 
angle occupied by such a hypothetical mid-latitude wind 
would preclude it from carrying a significant portion of the 
PNS's spin-down luminosity, which is primarily extracted 
by outflow originating near the equator (e.g., Bucciantini et 
al. 2006). 



4.4. Implications for GRBs 

We have shown that proto-magnetars born with P < P n w 
0.8 ms produce energetic, neutron-rich winds; it therefore ap- 
pears that the birth of very rapidly rotating proto-magnetars 
should produce neutron-rich GRBs. This is, however, not nec- 
essarily the case. Even strongly magnetized proto-magnetar 
winds are baryon-loaded at very early times following the 
launch of the SN shock because the PNS is very hot and 
its already substantial thermally-driven mass-loss is enhanced 
due to centrifugal slinging. In fact, a proto-magnetar requires 
several seconds to contract and cool to the point that a wind 
launched from its surface achieves the high magnetization a > 
10-100 required to explain the Lorentz factors inferred from 
GRBs. By the time a proto-magnetar cools to the point that its 
wind becomes ultra-relativistic, it may no longer rotate suffi- 
ciently rapidly to remain neutron-rich. 

In Figure [5] we explore this possibility quantitatively by 
showing the magnetization a of a PNS's outflow (eq. |fT31 ) as 
a function of the PNS rotation period P from evolutionary cal- 
culations of a cooling, spinning down proto-magnetar with an 
initial period Pq = 0.6 ms < P a and for three different fixed sur- 
face dipole magnetic field strengths: 10 B„ p = 10 15 G, 3 x 10 15 
G, and 10 16 G. These evolutionary calculations, which are de- 
scribed in more detail in MTQ07 §4.1, assume a PNS cool- 
ing evolution similar to that given in equation ([8]) and map the 
one-dimensional neutrino-heated monopole wind calculations 
of MTQ07 onto a more physical dipole geometry using the 
axisymmetric two-dimensional relativistic MHD calculations 
of Bucciantini et al. 2006 (hereafter B06). Although Figure[5] 
only shows calculations for one initial rotation period Pq = 0.6 
ms, the PNS's spin-down timescale is dominated by times at 
which it is slowly rotating; hence, the evolutionary tracks for 
PNSs born with Pq < 0.6 ms are similar to those shown in 
Figure [5] for P > 0.6 ms, following a brief initial spin-down 
phase. 

Figure |5] shows that PNSs born with P < P n and > 3 x 
10 15 G spin down to a relatively neutron-poor state (P > P n ) 
in just a few seconds, before cooling sufficiently to produce 
an ultra-relativistic (cr 3> 1) outflow; these proto-magnetars 
could therefore not produce neutron-rich GRBs. On the other 
hand, a PNS born with fij p = 10 15 G spins down more slowly 
and is thus capable of producing an ultra-relativistic, neutron- 
rich outflow ~ 10-30 seconds following core bounce. Of the 
proto-magnetar' s initial rotational energy of E Iot ~6x 10 52 

10 The minimum stable neutron star rotation period P mm is uncertain theo- 
retically because it depends on the uncertain supranuclear density equation of 
state (EOS); depending on the EOS, detailed studies find that P m ;„ = 0.53—1.7 
ms for M ~ 1.4M (Cook, Shapiro, & Teukolsky 1994, Table 8). Thus, al- 
though we take the specific value Pq = 0.6 ms in our calculations to illustrate 
the conditions under which proto-magnetars can produce neutron-rich GRBs, 
if Pmm > P n ~ 0.8 ms the conclusions of this section are much simpler: proto- 
magnetars cannot produce neutron-rich outflow. Conversely, the detection of 
a neutron-rich GRB outflow from a confirmed magnetar birth would provide 
a constraint on the EOS. 
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FIG. 5 . — Magnetization a of a PNS's outflow (eq. 1151 ) as a function of the 
PNS rotation period P from evolutionary spin-down calculations performed 
during the Kelvin-Helmholtz cooling of a proto-magnetar with initial rota- 
tion period Pq = 0.6 ms at a time t = 1 s after core bounce and for different 
fixed surface dipole magnetic field strengths: B„ p = 10 15 G, 3 X 10 15 G, and 
10 16 G; for reference, times t = 1, 3, 10, and 30 seconds after core bounce 
are shown along the evolutionary tracks. So long as the rotation period of 
the proto-magnetar remains less than P ~ P n pa 0.8 ms its outflow remains 
neutron-rich (see eq. 1161 ). The proto-magnetars with B„ p = 3 X 10 15 G and 
10 16 G spin-down to a relatively neutron-poor state in just a few seconds, 
before cooling sufficiently to produce a 3> 1 outflows. In contrast, the proto- 
magnetar with i?J p = 10 15 G spins down less rapidly, allowing the PNS to 
produce simultaneously neutron-rich and ultra-relativistic outflow ~ 10-30 
s following core bounce. 

ergs, a substantial portion (~2x 10 52 ergs) is extracted when 
a < 10-100; although this early wind is potentially too 
mass-loaded to produce a GRB, it may enhance the energy 
and alter the morphology of the accompanying SN shock 
(TCQ04; MTQ07; Bucciantini et al. 2007a). The remain- 
ing ~ 4 x 10 52 ergs, part of which could produce a neutron- 
rich GRB, emerges when a > 10-1 00 on a somewhat longer 
spin-down timescale tj, which is ~ 10 2 - 10 3 s, depending on 
the fraction of the proto-magnetar's surface threaded by open 
magnetic flux (B06). 

In $3] we emphasized the need for magnetic acceleration to 
produce neutron-rich outflows from compact objects. Despite 
the absence of magnetic fields, however, Dessart et al. (2006) 
find neutron-rich (Y" ^0.25), purely neutrino-driven outflows 
in 2D radiation-hydrodynamical simulations of the very early 
evolution (f < 800 ms after core bounce) of a rapidly rotating 
neutron star newly-formed following the accretion- induced 
collapse (AIC) of a white dwarf. Because these calculations 
show that very rapidly rotating PNSs can produce neutron- 
rich outflows without the aid of magnetic acceleration, it may 
appear that neutron-rich GRB outflows could be possible un- 
der less restrictive conditions than described in this section. 
The purely neutrino-driven neutron-rich outflows found by 
Dessart et al. (2006) are only possible, however, because at 
very early times following core bounce the PNS is inflated 
with respect to its final, cooled radius and because a PNS ro- 
tating near break-up is strongly oblately deformed. Because 
the equatorial neutrinosphere radius R v can therefore exceed 
~ 100 km at early times (see Dessart et al. 2006, Fig. 7), the 
gravitational binding energy of a nucleon on the PNS sur- 
face at a moderately low latitude is Eb ~ 10(7?„/100km) _1 
MeV, comparable to the mean energy of the neutrinos driving 



the outflow (typically w 10- 15 MeV during the PNS's early 
cooling phase). Thus, despite the fact that neutrino absorp- 
tions favor a proton- rich composition (Y" > 0.5) because 
the electron neutrino flux dominates the electron antineutrino 
flux during early deleptonization (especially in the presence 
of rapid rotation; Thompson, Quataert, & Burrows 2005), 
less than a single neutrino is required to unbind a typical nu- 
cleon (i.e., <fi D < 1; eq. fT8l ); matter driven from low latitudes 
can therefore partially retain the neutron-rich composition of 
the PNS surface. Although purely neutrino-driven neutron- 
rich outflows are thus possible from PNSs at very early times 
(< 1 s) after core bounce, as the PNS cools and its radius 
shrinks (f> n will increase; purely neutrino-driven neutron-rich 
outflows will therefore not be possible after a few seconds 
following core bounce, even if the PNS remains distorted by 
continuing to rotate near break-up. In particular, once the 
PNS has cooled and contracted sufficiently to produce a GRB 
outflow, neutron-rich winds will only be possible if they are 
magnetically-driven and satisfy the constraints described in 
this section. 

5. ACCRETION DISK WINDS 

In this section we describe calculations of the structure 
and neutron content of axisymmetric, one-dimensional MHD 
winds launched from NDAFs. As shown schematically in 
Figure [6] we perform our calculations along a spherical, 
monopole flux tube centered about the position "C" a distance 
Rotan9 directly below the black hole along the disk rotation 
axis, where Ro is the distance from the black hole to the wind's 
launching point just above the accretion disk midplane and 
9 £ [0° , 90°) is the angle that the flux tube makes with respect 
to the midplane. The distance from a given position "P" along 
the outflow to the black hole and monopole center are denoted 
r and s, respectively. The physical solid angle of the wind Af2 
is chosen such that the conserved mass outflow rate is given 
by M = Aflpv p s 2 = 47rpoVp,o/? 2 ) so that A57 = 47rcos 2 9, where 
v p is the poloidal wind velocity, p is the wind's density, and 
'0' denotes quantities evaluated at the base of the outflow; 11 
note that for equatorial outflow (9 = 0) M reduces to the def- 
inition used in S|4] for proto-magnetar winds. Although this 
choice for A57 is somewhat arbitrary, the quantities of most 
interest, a and Yf, do not depend on our normalization for M. 

As in the proto-magnetar case, we assume an open poloidal 
magnetic field B p = Bq(sq/s) 2 , where sq = Ro/ cos (9 and Bo is 
the strength of the poloidal field at the outflow's base. As 
discussed by Levinson (2006), whose formulation and geom- 
etry are similar to ours, a more consistent approach would 
be to seek a self-similar solution to the trans-field equation 
(e.g., Li et al. 1992; Contopoulos 1994). However, simulta- 
neously including the wind's slow point topology and conse- 
quent neutrino-heated mass-loss in such a formalism is diffi- 
cult without multi-dimensional neutrino-heated MHD calcu- 
lations, a formidable numerical challenge. In addition, the 
large-scale poloidal field threading NDAF disks is uncertain 
and so a more detailed calculation does not seem warranted. 

As discussed in the proto-magnetar case, Y° is generally de- 
termined fairly close to where the wind is launched near the 
disk midplane and so the precise radial scaling of B p (oc s~ 2 ) 

1 1 We use '0' to denote quantities at the base of the outflow instead of V 
(as was used in because, unlike in the PNS case, the disk midplane can 
be neutrino transparent and therefore may not posses a neutrinosphere. 
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FIG. 6. — Geometry of NDAF disk winds. The wind is calculated along 
a flux tube that is inclined at an angle 9 with respect to the accretion disk 
midplane. The calculation is started from a location just above the disk mid- 
plane a distance Rq from the central black hole and at a distance sq from the 
monopole center (position "C"). A generic point "P" along the outflow is lo- 
cated at a distance r from the central black hole and a distance s from the 
monopole center. 

is likely less important than its magnitude near the base of the 
wind ps Bq. The properties of the wind, however, can depend 
strongly on 9. In particular for 9 < 60° matter that co-rotates 
at the Keplerian orbital velocity which is perturbed from the 
disk along the field line is unstable and can escape to infin- 
ity, even without additional energy deposition (Blandford & 
Payne 1982). This results in significant mass-loading for field 
lines with 9 < 60°; in addition, the field develops a substan- 
tial base toroidal component B^q. For 9 < 60°, NDAF winds 
are therefore slow, heavily mass-loaded, have B^o 3> Bo, and 
possess properties that are relatively independent of 9. On the 
other hand, more vertical flux tubes (60° < 9 < 90°), which 
resemble equatorial proto-magnetar winds more than equa- 
torial disk winds, have lower M, accelerate to faster asymp- 
totic speeds, and possess properties that depend sensitively 
on 9; in particular, the wind becomes increasing relativistic as 
9 — > 90° (Daigne & Mochkovitch 2002), and ultra-relativistic 
outflows from nearly vertical field lines may be capable of 
producing GRBs. 

In terms of Z? and M the magnetization (eq. ifPTl ) of an 
NDAF wind launched from radius Rq is given by 

a = B 2 Xni/Mc 3 , (21) 

where we have assumed that the base of the outflow rotates at 
approximately the Keplerian rate 51 K = (GM/Rq) 1 ' 2 . The rate 
of angular momentum lost through a non-relativistic (er < 1) 
wind is given by J w = MR\fl K , while for relativistic (a > 1) 
winds Ra — Rl and Jw — £mag/^K — <tMc 2 /J7k, where E mlig 
is the wind's energy loss rate. The ratio of Jw to the disk's 
angular momentum loss rate Jo = MdRq^k is therefore given 
by 

J J- ~ (MR 2 /M D R 2 ) = (Ms 2 /M D s 2 ) : a < 1 
Jd 

~ (B 2 R 2 /M D c) :a>l (22) 
Although a wind with Jw < Jd is physical because stresses 



internal to the disk can transfer angular momentum outwards, 
allowing accretion at Mo to proceed, a wind for which Jw > 
Jd is unphysical because it violates conservation of angular 
momentum. 

5.1. Numerical Procedure 

We calculate NDAF wind solutions starting from a location 
just above the accretion disk midplane at a fixed cylindrical 
radius Rq = 14-R g w6x 10 6 cm, which is roughly the loca- 
tion of the disk's peak neutrino flux. The wind's local neu- 
trino luminosities and energies are determined from CB07's 
marginally neutrino optically thin Md = 0.2M Q s" 1 , a = 0.03 
NDAF solution. We focus on winds driven from low M D , 
neutrino transparent disks because low M D accretion, due to 
its lower neutrino-driven mass-loss, is probably more likely 
to produce neutron-rich, ultra-relativistic outflow. When con- 
sidering NDAFs with M D < 0.2M Q s" 1 (down to M D ~ M ign ; 
et l- CQX we scale the boundary conditions and neutrino emis- 
sion properties analytically from the Md = 0.2M Q s" 1 solu- 
tion, as discussed below. 

We use Newtonian gravity, including just the point mass of 
a M = 3Mq black hole. The evolution equations for v p , v^, 
Z?0, p, and T are the same as those used in the PNS case with 
r replaced by s, except that the gravitational acceleration in 
the radial momentum equation (MTQ07, eq. [7]) is now the 
component projected along the field line: 

= GM(s -s sm 2 9) 

8s (s 2 + sm 2 9[s 2 -2 SSo ])V 2 ' 

Our microphysics is essentially identical to that in the PNS 
case. In particular, the electron fraction is evolved according 
to equation (fl4] i, again with r replaced by s; note that because 
of the disk's significantly lower density, degeneracy effects 
are less important near the base of the wind than in the PNS 
case. The electron and anti-electron neutrino energy fluxes 
(F Uii ,Fd c ) used to calculate the neutrino absorption and heating 
rates near Rq are computed from CB07's a-disk solution in 
the same manner as the local neutrino energy densities that 
were used to calculate Y£ for Figure [T] to within a factor of 
~ 2 the total spherically-equivalent neutrino luminosity obeys 
L v = L Ve +L 9e ee 4%R 2 (F Ke +F„ e ) « 7yM D c 2 , where 77 « 0.04 
is the total radiative efficiency for a non-rotating black hole 
from CB07. In a similar manner as for PNSs, we take the 
neutrino flux intercepted by the outflow to be approximately 
constant with radius for r < Ro and to decrease cx r~ 2 for r > 
R (T01 eq. [24]; see Surman & McLaughlin 2004, 2005 for 
a more detailed treatment). Although the geometry that we 
assume is simplistic, Y" is set relatively close to the base of 
the wind and thus depends primarily on the neutrino energy 
density in the vicinity of ^o- Indeed, as Figure [3] illustrates, 
the advection and pair capture rates rise and fall above the 
base of the outflow, respectively, much more rapidly than the 
neutrino absorption rates decrease; this would likely remain 
true for any realistic geometry. 

Unlike PNSs, NDAFs do not produce significant tau or 
muon neutrino emission and so annihilations and neutral- 
current interactions from these neutrino species are ignored. 
The mean electron and anti-electron neutrino energies near 
Rq, which are not expected to vary strongly with Md for neu- 
trino optically thin accretion (see eq. |9l), are taken from Fig- 
ure[l]as (e„„) = 10.5 MeV and (e Pt ) = 13.1 MeV; fromCB07's 
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solution we also determine that Fp e ~ 1 .2F Ve near Rq so that 
Y" w 0.51 (see eq. (6) and Fig. [T}. We set the wind's inner 
density po equal to the disk midplane density po if the elec- 
tron neutrino optical depth along s to the midplane is r„ < |; 
otherwise, we instead choose /5„ to enforce t„ ~ | at so- 

Although viscous heating generally dominates neutrino 
heating in the midplane of NDAFs, recent radiation MHD 
simulations suggest that little energy is dissipated in the disk 
corona (Turner 2004; Hirose, Krolick, & Stone 2006; Blaes 
et al. 2006b; Krolik et al. 2006). For this reason, we neglect 
viscous heating in the wind entirely and set the base temper- 
ature To of the wind by balancing neutrino cooling with just 
neutrino heating, as in the PNS case; note, however, that sig- 
nificant viscous heating in the wind would likely result in both 
additional mass-loss and deneutronization. The disk's mid- 
plane temperature T D , which is set by the balance between 
viscous heating and neutrino cooling (see eqs. ll9l- fT0l ), is 
therefore generally higher than 7b. 

Although NDAFs are efficiently cooled and geometrically- 
thin, radial pressure support is not completely negligible and 
so the disk's angular rotation frequency fi, which we use to set 
the wind's inner angular velocity, is slightly sub-Keplerian: 
il 2 = n^l-H 2 /^), where H w 0.2R is the disk scale height 
near Rq, which, like To, is approximately independent of Mo 
for neutrino optically thin accretion. 

Our disk wind calculations can be compared with sim- 
ilar one-dimensional flux tube calculations by Daigne & 
Mochkovitch 2002 (DM02), Pruet et al. 2004 (P04), and 
Levinson 2006 (L06). DM02 calculated the requisite condi- 
tions for ultra-relativistic outflow from hyper-accreting disks, 
including neutrino heating and cooling and, in the neutrino 
optically thin case, a simplified viscous heating prescription. 
Because DM02 was primarily concerned with obtaining M as 
a function of disk conditions, they only considered wind con- 
ditions near the sonic point; in addition, DM02 assumed co- 
rotation rather than accounting for the magnetic field explic- 
itly. L06 improved upon the calculations of DM02 by includ- 
ing the full equations of general relativistic MHD. Although 
L06 explored the effects of finite Bo, L06 concentrated, like 
DM02, on the sub-slow magnetosonic regime and his calcu- 
lations did not capture the Alfven or fast magnetosonic radii. 
Although this approach allowed L06 to calculate M as a func- 
tion of the open magnetic flux and L v , the base toroidal field 
fi^.o remained a free parameter in L06's formulation. Because 
B^fi is associated with the conserved magnetic induction it 
is fixed in our calculations by the fact that our steady-state 
winds pass smoothly through all three MHD critical points. 
L06 speculated on the potential deneutronization of magne- 
tized NDAF winds by noting the similarity between the ad- 
vection and relevant weak interaction rates; he did not, how- 
ever, calculate the evolution of Y e explicitly. 

Finally, P04 investigated nucleosynthesis in collapsar disk 
winds by solving the equations of hydrodynamics and by 
evolving Y e from the disk midplane. P04 included neu- 
trino heating and viscous heating through an a-prescription; 
although they neglected neutrino absorptions in evolving 
Y e , they also argued for the generic deneutronization of 
thermally-driven winds. Although P04 did not include the 
effects of magnetic fields on the wind explicitly, their out- 
flows were calculated along well-defined, vertically-directed 
flux tubes and were artificially forced to co-rotate outside the 



base of the wind, presumably to mimic the effect of a strong 
poloidal magnetic field. Our calculations are, to the best of 
our knowledge, the first to fully calculate the effects of MHD 
on the evolution of Y e in NDAF outflows and the first to cap- 
ture all three MHD critical points. 

5.2. Numerical Results 

The solid lines in Figure UJ show calculations of the mass- 
loss rate M, asymptotic electron fraction Y", magnetization a 
(eq. ETl ). and angular momentum loss rate Jw compared to 
that required for accretion Jo (eq. l22l ) for NDAF winds with 
base poloidal field B Q = 10 13 G, 10 14 G, and 10 15 G as a func- 
tion of flux tube angle 9 for the M D = 0.2 M© s" 1 , a = 0.03, 
M = 3M Q , a = NDAF solution of CB07 (the solution on 
which Figure[T]was based). Because magnetic pressure at the 
base of the outflow exceeds the total thermal pressure in the 
midplane of the disk for Bo>2x 10 15 G, fields much stronger 
than those we have considered in Figure UJ are probably un- 
physical because they could not be self-consistently anchored 
to the disk. Our disk wind calculations are summarized in 
Table |3] 

Figure UJ shows that M is very large and relatively inde- 
pendent of 9 for 9 < 60° but that for 9 > 60° M decreases 
rapidly with increasing 9. The mass-loss rate also increases 
with increasing Bo, except for the largest angles, for which M 
saturates for sufficiently large Bo, no longer increasing with 
increasing Bo (as in the PNS case; see Fig. |4). There is no 
similar saturation for smaller 9 because the large values of M 
preclude the outflow from co-rotating out to the sonic point 
for physical Bo- 12 

Figure UJ also shows that although Y" is relatively low (~ 
0.3) for the high-Bo, low-0 wind solutions, all of the solutions 
in FigureUJwith Y" <C Y" ~ 0.51 are unphysical because they 
have Jw > Jd and thus extract angular momentum at a rate 
exceeding that required for accretion through the disk from 
which the wind's boundary conditions were derived; in fact, 
invoking the criterion Jw < Jd (see eq. 1221 ). the only physical 
solutions in Figure UJ are those with Bo= 10 13 G for all 9 and 
Bo = 10 14 G for 9 > 80°. Also note from Figure UJ that, in 
addition to having Y" ~ 0.5, all of the physical solutions for 
Mo = 0.2Mq s" 1 are non-relativistic (er < 1). Because our 
calculations have thoroughly spanned the physical parameter 
space of Bo and 9, this conclusion is robust, at least for Mo = 
O.2M s _1 . 

For lower Mo the deneutronizing neutrino luminosity and 
neutrino-driven mass-loss rate are lower and thus relativistic, 
neutron-rich outflow may be more likely. To explore this pos- 
sibility we have also calculated wind solutions for outflows 
driven from an NDAF with M D = 10" 2 M Q s" 1 ~ M iga (a = 
0.03), keeping the other parameters identical to those of the 
higher M D case (a = 0.03, M = 3M Q , R = l4R g ). To compute 
these low Mo solutions we decreased the neutrino fluxes by 
a factor of 20 from our Mo = O.2M s" 1 calculation, left the 
mean neutrino energies unchanged (they are not expected to 

12 Because the base of the wind rotates at a slightly sub-Keplerian rate 
due to radial pressure support in the disk, even for 8 < 60° mass-loss would 
saturate atMsj 47ri? 2 l poCj i o ~ 10- 100M© s~' for sufficiently large Bo (~ 
10 16 G; the co-rotating limit), where c s o ~ 0. 1 c is the sound speed near the 
base of the outflow. 



13 



10 



10 1 ■ 



io- 2 ■ 



10 : 



■ 






■ +■■ 








10 13 G 




+ 


10 u G 




. o 


10 15 G 






M D = 0.20 M s" 1 






M D = 0.01 M s 1 

















50 



60 70 
8 (degrees) 



80 



90 



0.5 r 



0.4 r 



> 0.3 : 



0.2 r 



0.1 

























— *~Z?y 






































+ 


























4 























50 



60 70 
(degrees) 



80 



90 




60 70 
9 (degrees 



60 70 
9 (degrees) 



90 



FIG. 7. — Mass loss rate M, asymptotic electron fraction F", magnetization a (eq. 1211 ). and the ratio of the wind to accretion angular momentum loss rates 
Jw/Jd (eq. 1221 ) for NDAF winds as a function of the angle 9 between the wind's flux tube and the disk midplane (see Fig. [6}; base poloidal magnetic field 
strengths B = 10 13 G (asterisk), 10 14 G (cross), and 10 15 G (diamond) are considered. The solutions were calculated starting from a location just above the 
disk midplane at a radius So = 14fi g ~ 6 X 10 6 cm from the central M = 3Mq, non-rotating (a = 0) black hole and employed boundary conditions derived from 
the a— disk NDAF solution of Chen & Beloborodov (2007) with a = 0.03 and accretion rate Mo = 0.2Mq s (solid lines). Also shown are wind solutions 
for Md ~ 1O~ 2 M0 s and Bo = 10 14 G. For 9 < 60° the wind properties are relatively independent of 9 while for 9 > 60°, M decreases and a increases 
rapidly with 9. For the most vertical field lines (9 — * 90°), along which the outflow experiences minimal centrifugal support, M approaches its purely thermal, 
neutrino-driven value (eq. 1241 ). Of the solutions shown, only those with Jw/Jd < 1 are physical. All of the physical solutions for Mo = 0.2Mq s _1 are 
non-relativistic (a < 1) and do not possess an asymptotic neutron excess (Y" fa ~ 0.51). Winds from the Mo = 10~ 2 Mq s~' NDAF have lower M and larger 
a than the winds driven from the higher Mo disk; however, the solutions for Mo = 10~ 2 Mq s~' with Jw < Jd still have Y" ~ 0.5. 



vary strongly with Mn; see eq. J5)), and decreased the disk 
midplane density po ~ po oc Md- We only calculated solu- 
tions with 5o^3x 10 14 G because for larger field strengths 
the magnetic pressure would exceed the thermal pressure in 
the disk. 

In Figure [7] we show the Mo = 1O~ 2 M0 s" 1 calculations 
for Bo = 10 14 G with a dotted line for comparison with the 
higher Md solutions. Although the properties of these low 
Md winds and their dependence on 9 qualitatively resemble 
the higher Md solutions, Y" and M are generally lower and a 
is higher than for the winds driven from the more neutrino- 
luminous disk. Most notably, the high-0 solutions are now 
both physical (Jw < Jd) and relativistic (a > 1), illustrating 
that nearly vertically-directed outflows from low-Mo NDAFs 
represent plausible GRB central engines. However, these 
high-^, relativistic winds still possess no significant neutron 
excess (Y" ~ 0.5); this means that, even for Md ~ Mi Sn , si- 



multaneously neutron-rich and ultra-relativistic outflow ap- 
pears unlikely. One might think that Md (and hence L v ) could 
be lowered further until neutron-rich, ultra-relativistic outflow 
was possible at some intermediate latitude; a key difference, 
however, between the NDAF problem considered here and 
the PNS problem considered in f|4]is that, unlike in the PNS 
case, the disk's neutrino luminosity cannot be decreased be- 
low L uhn = r)M ign (a = 0.03)c 2 - 10 51 ergs s" 1 (eq. Q). For 
Md < M- isn the disk is geometrically thick, and the mid plane 
is not necessarily degenerate and neutron-rich (see ^5. 5b . 

5.3. Analytic Constraints 

As in the PNS case, the purely thermal, neutrino-driven 
mass-loss rate from NDAFs can be estimated analytically by 
requiring that the energy used to unbind a typical nucleon Eb 
be supplied entirely by neutrino heating; equation (|T9b can 
therefore also be used to estimate M± for NDAF disk winds, 
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provided that /?„ — > Rq, (e„) is taken from equations (O and 
( [Tol l, and that, as discussed below, we include the effects of 
the lower effective gravity g e s near the base of the outflow 
than in the PNS case. For neutrino optically thin accretion we 
find that M th for NDAFs is given by 



1(T 



Mian 



5/3 



f gzs/iGM/Rl) 



H/R 



,.119/90 301/90 

M 3 a 0A < 



"p 
6R e 



6R„ 



79/60 



M© s 



(24) 



where we have used equation (JTJ to scale the neutrino lumi- 
nosity to the value L„,i gn associated with the ignition accretion 
rate M ign (eq. fl]) because NDAFs only exist for M D > M ign . 

The purely neutrino-driven mass-loss rate for PNSs 
(eq. fT9l ) was derived for outflows emerging antiparallel to 
the PNS's gravitational field and by assuming substantially 
sub-break-up rotation. NDAF winds, however, are driven 
from the disk midplane at an angle inclined with respect 
to the black hole's gravitational field (i.e., 9^0) and even 
non-magnetized centrifugal support can be important near the 
base of the outflow because the disk rotates at nearly the 
Keplerian rate. This means that the effective surface grav- 
ity g e ff(0), which is the difference between the gravitational 
acceleration inward along s {g s (9); eq. l23l ) and the cen- 
trifugal acceleration outward along s (vi/i), is always less 
than the gravitational acceleration for purely equatorial out- 
flow (g s (6 = 0) = GM/Rq) that was used to derive M m for 
PNSs. Assuming that the wind conserves angular momen- 
tum away from sq (i.e., = HR^(s / sq)~ 1 ) and concentrating 
on the wind's inner, quasi-hydrostatic atmosphere, we find 
that g e f t /(GM /R^) is roughly constant with 9 and is approxi- 
mately equal to (H/R ) w 0.3a l) / ^M~ l/l0 (R /6R g )y 20 , where 
H is the NDAF's vertical scale height near/?o; this calculation 
justifies our characteristic scaling for g e ff in equation (124) . 

Although equation d24l i represents the thermal, neutrino- 
driven mass-loss rate from an NDAF, only the nearly-vertical, 
9 = 89° wind solutions shown in Figure [7] (which experience 
minimal centrifugal support) in fact reach M ~ Mth- Thus, 
equation (124-b should be taken as the absolute minimum NDAF 
mass-loss rate; in the presence of even a modest magnetic 
field with a non-vertical inclination, M is significantly larger 
that M m . Because NDAFs possess an absolute minimum 
mass-loss rate, their outflows also posses an absolute maxi- 
mum magnetization cr max , which is given by 
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where (3 P is the ratio of the disk's midplane thermal pressure 
Pd to the magnetic pressure associated with the wind's base 
poloidal field B\/%n and (3^, is similarly defined for the disk's 
midplane toroidal magnetic field B^o- In equation d25l l we 
have written a mwi in terms of (a/^) -1 because for angular mo- 
mentum transport via the MRI, local shearing box simulations 
find that a(3^, ~ 1 (Hawley, Gammie, & Balbus 1995); we also 



write equation (l25l l in terms of P<p/(3p = B\jB 1 bD because if 
the open poloidal field is generated through a dynamo from 
the toroidal field this ratio is unlikely to exceed unity (and is 
probably much less). 

Equation ((25) shows that cr max oc M D 2 ^ , which explains why 
the low Mp solutions in Figure [7] produced more relativistic 
outflows than the high Mo solutions under the physical con- 
straint (3 P > 1 that we imposed on our calculations. Equation 
( |25] > also shows that NDAFs can, in principle, produce ultra- 
relativistic (ct 3> 10 2 - 10 3 ) outflows from small radii (Rq ~ 
^isco), provided that a significant fraction of the magnetic 
energy present in the disk is associated with a large-scale, 
open poloidal field (i.e., /3 p ~ (3 t ). For a = 0.1, M = 3M Q , 
/5 P = (3$, ctftcf, = 1, and for accretion onto a rapidly rotating 



black hole (so that Rq «S p R! 2R g ) equation (125 D gives a 
10 2 (M D /M ign )- 2 / 3 . In this case ultra-relativistic outflow could 
occur for Mo ~ M gn and could thus accompany a substantial 
accretion power E acc = r)M D c 2 ~ r]Mi gn (a = 0. l)c 2 ~ 10 51 ergs 
s _1 . Because £,cc represents the maximum MHD luminosity 
of a GRB-producing jet, this shows that, under ideal condi- 
tions, disk winds from NDAFs form plausible GRB central 
engines. 

Equation (f25j) shows that under some circumstances ultra- 
relativistic outflow from NDAF disks is plausible. How- 
ever, the maximum magnetization for material with a neutron 
excess is lower than er max because the minimum mass- 
loss rate that must accompany a neutron excess is substan- 
tially larger than M t h according to the same arguments that 
were used in ^4.31 for the PNS case; indeed, using equation 
d25l l and the NDAF analog of equation ( fT8l we find that 
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(26) 



Equation d26l > shows that neutron-rich, relativistic outflow is 
very unlikely from NDAF disks. For Mo = M gn , a = 0.1, 
M = 3M Q , P p = f3<p, a/?^ = 1, and for outflow launched from 
near /?j sco of a rapidly spinning black hole (Rq = R p = 2R g ), 
equation d26i > gives <r max ~ 10, insufficient to explain the ultra- 
relativistic outflows inferred from GRBs. 

If the fiducial scalings of equation d26j > are not adopted, 
ultra-relativistic outflow with a neutron excess may be pos- 
sible in some circumstances. For instance, if a = 0.01 instead 
of a = 0.1 (again for M D = M gn , M= 3M Q , /3 p = (5$, a/3 ct> = 1, 
^o = = 2R g ) then cr max ~ 10 2 — 10 3 . In this case, however, 
£ acc (a = 0.01) ~ 10 5l) (M D /Mi gn ) ergs s _1 and so a very large 
fraction of the accretion energy would need to be deposited 
in the ultra-relativistic outflow to explain the observed lumi- 
nosities of GRBs. Likewise, although a(3 v <C 1 is possible if 
the poloidal field threading the disk is not the result of a local 
dynamo (instead resulting from, e.g., magnetic flux advected 
from large radii in the disk; e.g., Spruit & Uzdensky 2005), 
the constraint Jw < Jd (eq- l22l ) independently requires that 

a/3 p > 2/3( J R„/^)- 1 / 2 (^o/H) « 4a" 1 / 'Vf ^,)" 17 / 20 ; 
thus, outflow launched from near R KCO must have [3 p a > 1 
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because otherwise it would carry away angular momentum at 
a rate exceeding that required for matter to accrete at the rate 
Mb- Lastly, we reiterate that the true thermally-driven mass- 
loss rate can, in principle, far exceed the purely neutrino- 
driven value M t h (eq. 1241 ) if viscous heating is important 
near the base of the wind; if this were the case, (rJJ lax would 
be substantially reduced below the value given by equation 
d26| i. We thus conclude that neutron-rich GRB outflows are 
unlikely from NDAF disk winds. 

5.4. Cross-Field Neutron Diffusion 

In the previous section we have argued that highly relativis- 
ti c C" 7 ^ 100- 1000) winds driven from the innermost radii 
of NDAFs are unlikely to be intrinsically neutron-rich. How- 
ever, free neutrons are uncharged and may collisionally dif- 
fuse across magnetic field lines to the polar region (hereafter, 
the "jet") from an adjacent, more baryon-rich wind (Eich- 
ler & Levinson 1999; Levinson & Eichler 2003, hereafter 
LE03; McKinney 2005b). If the total neutron mass diffu- 
sion rate MjJ lff dominates the mass-loading of the polar jet then 
the highly relativistic polar outflow will be significantly "pol- 
luted" by neutrons and may end up neutron-rich after all. 

Neutron diffusion into the polar jet from the adjacent mass- 
loaded wind is limited to a surface area ~ Ans^ip, where s a is 
the distance from the base of the wind to where free nucleons 
recombine into a-particles (which are charged and therefore 
cannot efficiently diffuse across field lines) and we have as- 
sumed that s a S> Ro/<p, where ip = tt/2-9 <§; 1 is the opening 
angle of the jet. Neutron diffusion is limited by elastic proton 
collisions, with a rate (<7 n _ p v re i) — 10~ 15 cm 3 s" 1 and a corre- 
sponding mean free path A n - P — Vth/(«p(cn-pV re i)), where n p 
and v t h ~ (kT/nin) 1 / 2 are the proton number density and the 
ion thermal speed, respectively. Following LE03, we assume 
that the density gradient length scale separating the mass- 
loaded wind and the axial jet is given by I ~ (vthTdynAn-p) 1 / 2 , 
where Td yn = s/v p is the wind's dynamical timescale. 13 Us- 
ing the Fick-diffusion number flux of neutrons into the jet of 
F n n n v t h(A n _ p //), we estimate that 

.. itT , . { Airs a kT a M\ 1/2 

M n dlff ~ 47r^> n F n L = , " 

\ (CTn-pVrel/ / 

VMeVy V10 8 cm/ \0AMqS~ 1 J 

(27) 

where M w 4irip 2 m n v p (n p + n n )s 2 , n„, and T a = T(s a ) are the 
wind's mass-loss rate, neutron number density, and temper- 
ature at a-particle recombination, respectively, and we have 
assumed that Y e ~ 0.5 in the wind. We evaluate equation d27l > 
at s a because MjJ lff is dominated by the largest radii at which 
the wind is still primarily free nucleons. 

For the relatively moderate entropies (S a < 10 2 
kg baryon" 1 ) that characterize neutrino-heated, 
magnetocentrifugally-driven winds (see Tables [2] and [3), 

13 This choice for / is appropriate for a very abrupt transition in the wind's 
density with cylindrical radius, such as between field lines threading the disk 
and those threading the black hole's event horizon. A perhaps more natural 
(but less conservative) choice for / is the cylindrical radius in the wind at 
a— particle recombination (?s <ps a ), which would produce an even smaller 
M* ff than is given in equation (27). 



a-particles form at a high temperature (T a ~ 1 MeV), which 
obtains relatively close to the base of the wind (s a < 10 7 
cm). In this case, even for M ~ M Q s" 1 , equation d27i > 
gives M„ lff < 10" 8 M o s" 1 . If the axial jet is itself driven 
from the disk, MjJ lff is thus significantly lower than the 
minimum mass-loss already supplied by neutrino heating 
(eq. 11241 ); hence, cross-field neutron diffusion is ineffective 
at segregating neutrons in low entropy NDAF winds. 

If, on the other hand, the axial jet is powered by v - v an- 
nihilation or the Blandford-Znajek process and has little or 
no intrinsic baryon-loading (such as if it threads the black 
hole's event horizon), then M* ff , although small, may dom- 
inate the jet's mass loading. For instance, for a polar jet 
power of E ~ 10 50 - 10 51 erg s" 1 equation ( |27T i shows that 
diffusive neutron mass-loading from an encasing wind with 
a mass-loss rate M ~ 10~ 2 M Q s" 1 would, by itself, limit the 
jet's asymptotic Lorentz factor to T ~ 10 4 - 10 5 . Asymp- 
totically neutron-rich outflow may result in this case if the 
jet remains "clean" to large radii; elucidating the observable 
consequences of such very high-r neutron-rich outflows will, 
however, require additional work. Lastly, we note that al- 
though our calculations show that NDAF winds probably pos- 
sess moderate entropy, previous works that have considered 
diffusion into the jet have focused on very high entropy out- 
flows characteristic of hydrodynamic "fireballs" in the GRB 
literature (LE03; McKinney 2005b). These calculations find 
larger Mj| lff than we have estimated in equation (|27T i in large 
part because a-particles do not form until much larger radii 
in high entropy winds. Furthermore, if Y" < 0.5 in the encas- 
ing baryon-rich wind and some neutrons remain free to radii 
larger than s a , may be larger than predicted by equa- 
tion d27b . Evaluating the free neutron fraction that remains 
in magnetically-driven disk winds following the freeze-out of 
any a-process or r-process capture that may occur following 
a-particle recombination is, however, beyond the scope of 
this paper. 

5.5. Thick Accretion Disk Winds 

For both proto-magnetars, where L v oc t~ l starting ~ 1 s 
following core bounce (eq. [8]), and NDAFs, where L v oc 
Mo oc f " 5 / 3 is expected at late times from "fall-back" accretion 
(Chevalier 1989; Woosley & Weaver 1995) or M D oc r a with 
a ~ 1 due to the viscous evolution of a disk of finite mass, the 
low neutrino luminosities attained at late times are the most 
favorable for ultra-relativistic, neutron-rich outflows. How- 
ever, while proto-magnetars must maintain sub-millisecond 
rotation for most of the Kelvin-Helmholtz cooling epoch to 
produce late-time neutron-rich GRB outflows, for NDAFs, 
which always rotate at nearly the Keplerian rate, the difficulty 
is more fundamental; for L v < L„.i gn (eq. Q) the disk is no 
longer efficiently cooled and the disk midplane may no longer 
remain dense, degenerate, and neutron-rich. For both collap- 
sars and the accretion accompanying compact object mergers 
there will thus come a time f; gn after which Ma < M; gn and the 
disk will transition from an NDAF to an advection-dominated 
thick disk (Narayan & Yi 1994; Narayan et al. 2001). 

Although NDAFs enter /3-equilibrium on an accretion 
timescale (B03a), thick disks generally do not. Using equa- 
tion (48) from B03a we find that /^-equilibrium is only es- 
tablished in a thick disk at radii smaller than a critical f3- 
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equilibrium radius Rp, which is given by 

Rp = 34(14)(M D /M ign ) 10 / 13 a^ /3 X /39 ^ (28) 
for accretion onto a black hole with spin a = 0(0.95). Because 
Rp < fiisco for M D < 0.1(0.03)M ign , equation (J28J shows that 
matter accreting through a thick disk at a rate Mo <C M{„ n is 
not in /3-equilibrium at any radius. The neutron content of 
the disk at small radii for t fi gn (and thus of any potential 
late-time GRB-producing outflow) will therefore depend on 
the composition of the matter feeding it. In particular, late- 
time GRB outflows from collapsar disks, which are continu- 
ally fed from large radii by their progenitor's stellar He enve- 
lope (which has Y e ~ 0.5), will not be neutron-rich. 

On the other hand, accretion disks formed from compact 
object mergers, which are usually fed from the tidal disrup- 
tion of at least one neutron star 14 (e.g., Rosswog & Ramirez- 
Ruiz 2002), are initially neutron-rich. However, the disks 
formed from compact object mergers are expected to be 
more compact than collapsar disks (with circularization radii 
< 10-30r? ? ) and are probably not continually supplied with 
substantial mass from large radii; the late-time neutron con- 
tent of thick disks from compact object mergers therefore 
depends on the evolution of Yf immediately following the 
NDAF to thick disk transition at fi gn . Because equation d28l 
shows that a thick disk with Mq < M- ign does have sufficient 
time to enter /3-equilibrium, the late-time electron fraction 
in a thick disk from a compact object merger depends on 
whether /3-equilibrium near R lsco in a thick disk for M ~ M; gn 
favors a neutron-rich or a proton-rich composition. Although 
NDAFs are always sufficiently dense and degenerate to fa- 
vor a neutron-rich composition, /3-equilibrium in a thick disk 
only favors n/p > 1 for radii smaller than a critical neutron- 
rich radius R n ; using B03a equation (50) we find that 

R n = 24(2)(M D /M ig „) 2 «o^3 4/3 ^ (29) 

for accretion onto a black hole with spin a = 0(0.95). Equa- 
tion (1291 shows that when a disk transitions from an NDAF 
to a thick disk, the matter near Ri sco may be driven to either 
a neutron-rich or a proton-rich state, depending on a, M, and 
the extent of the disk. Thus, although disk winds from NDAFs 
are unlikely to produce neutron-rich GRB outflows, neutron- 
rich outflows may be possible from the thick disks associated 
with compact object mergers at late times (f > fj gn ) or from 
collapsar disks at t ~ l lgn (i.e., after the thick disk transition 
but before the disk is fed by additional, neutron-poor material 
from large radii). 

6. DISCUSSION 

By calculating the structure and neutron content of 
neutrino-heated MHD winds driven from the neutron-rich sur- 
faces of proto-magnetars and NDAFs, we have delineated the 
conditions under which a large neutron excess can be pre- 
served in these outflows. We have focused on the conditions 
for simultaneously neutron-rich and ultra-relativistic outflows 
because magnetized winds from hyper-accreting disks and 
newly-formed magnetars are plausible GRB central engines; 
despite being difficult to distinguish on other grounds, each 
of these central engines may possess a distinctive nucleonic 

14 An exception are the mergers of black hole-white dwarf binaries, which 
may produce long-duration GRBs (Fryer et al. 1999b). 



signature. If the consequences of neutron-rich GRB outflows 
enumerated in SfTJcan be identified or constrained, magnetar 
and black hole models for GRBs may thus be observationally 
distinguishable. 

Although GRB central engines are often neutron-rich (Pruet 
et al. 2003; B03a), we find that ultra-relativistic neutron-rich 
outflows are possible only under surprisingly limited circum- 
stances. Central engines that are sufficiently dense and degen- 
erate to be neutron-rich must be efficiently neutrino-cooled. 
For the resulting sub-virial temperatures, several of the ther- 
mal neutrinos released by the central engine must be absorbed 
by a typical nucleon for it to escape the deep gravitational po- 
tential due to neutrino-driving alone. Since neutrino absorp- 
tions from efficiently neutrino-cooled central engines usu- 
ally favor an asymptotic electron fraction Y" > 0.5, purely 
neutrino-driven outflows are generally driven back to a rel- 
atively neutron-poor state, with the precise value of Y" de- 
termined by the neutrino spectrum of the central source (see 
eq. ]6) and surrounding discussion). 

Additional forces (e.g., magnetocentrifugal) can prevent 
deneutronization by supplying most of the binding energy 
needed to escape the central engine's gravitational potential 
well. However, Y° is set so close to the base of the outflow 
that the very inner, hydrostatic atmosphere of the wind must 
be altered for Y" <C 0.5 to obtain; this unavoidably increases 
the wind's mass-loss rate M (see eq. lfT8l ). Indeed, a generic 
anti-correlation between M and Y° is evident in our numerical 
calculations shown in Figures |2] Hi and [7] Neutron-rich GRB 
outflows are thus difficult to produce because the minimum 
neutron-rich mass-loss rate often precludes ultra-relativistic 
speeds. Other heating (e.g., viscous) that may be present in 
the outflow in addition to neutrinos (which must be present) 
will only further increase M and, through additional entropy 
deposition and pair creation, further deneutronize the outflow. 

Our conclusion that simultaneously ultra-relativistic and 
neutron-rich outflow is difficult to produce depends on the as- 
sumption that Y" is not significantly less than 0.5. Therefore, 
to be more precise, the conclusion of this paper is that the nu- 
cleonic content of ultra-relativistic outflows driven from effi- 
ciently neutrino-cooled central engines is typically set by an 
equilibrium with neutrino absorptions. Thus, if the neutrino 
spectra and luminosities of PNSs and NDAFs are significantly 
different from what current calculations find, and Y% <C 0.5, 
neutron-rich outflows from GRB central engines may be more 
common. The next step in improving our understanding of 
magnetocentrifugal winds from GRB central engines is to in- 
clude the effects of a strong magnetic field on the neutrino 
interactions and the equation of state for the leptons. For the 
parameters considered in this work, the latter should be more 
important, modifying the pressure and entropy profiles. 

With these general constraints in mind, we now discuss the 
prospects for neutron-rich outflows from individual central 
engines. Our conclusions are summarized in Table [4] 

6.1. Proto-Magnetars 

Proto-magnetars with surface magnetic field strengths B v > 
10 14 — 10 15 G can produce neutron-rich outflows, but only 
for rotation periods P < P n as 0.8 ms (see Fig. 0. If the 
minimum stable neutron star rotation period P m ; n exceeds 
P n , then neutron star birth should not be accompanied by 
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substantially neutron-rich outflow. 15 On the other hand, if 
P n < P m ; n and steady-state neutron-rich winds from magne- 
tar birth are indeed possible, they should be restricted to 
events with a total GRB plus SN energy exceeding E mt (P n ) w 
4 x 10 52 ergs (absent significant gravitational wave losses). 
Such extremely energetic events are likely rare, even among 
GRBs. Furthermore, not all magnetar births that release 
> 4 x 10 52 ergs will necessarily produce simultaneously ultra- 
relativistic and neutron-rich outflows because proto-magnetar 
winds are heavily mass-loaded at early times following core 
bounce. If the proto-magnetar's dipole field strength exceeds 
B^ p ~ 10 15 G, spin-down is so rapid that neutron-rich out- 
flow is unlikely by the time the PNS has cooled sufficiently 
that a > 100- 1000 (see Fig. |5). Since the PNS spin-down 
power for P w P a and B^ p w 10 15 G is - 3 x 10 49 ergs 
s" 1 (Spitkovsky 2006; B06), neutron-rich GRBs from proto- 
magnetars possess a maximum beaming-corrected luminosity 
L 7 ~ 10 49 (e/0.3) ergs s , where e is the efficiency for con- 
verting outflow to gamma-ray energy. 

6.2. NDAFs 

Although neutron-rich GRB outflows are possible from a 
subset of proto-magnetars, we find that they are unlikely to 
originate from NDAFs under any circumstances. In agree- 
ment with previous works (e.g., Daigne & Mochkovitch 2002; 
Levinson 2006), we find that outflows with a > 100 are pos- 
sible from the innermost radii of NDAFs around rapidly rotat- 
ing black holes (a « 1), despite the minimum neutrino-driven 
mass-loss rate (eq. 11241 ). The significantly larger M required 
for a large neutron excess, however, precludes neutron-rich 
NDAF outflows from attaining cr > 1 — 10 for significant ac- 
cretion powers (eq. l26l ). Furthermore, because modest en- 
tropy, magnetocentrifugally-driven winds only possess free 
nuclei over a relatively limited range of radii, cross-field dif- 
fusion is ineffective at polluting otherwise neutron-poor axial 
jets with free neutrons from adjacent, more heavily baryon- 
loaded winds (see $5.41 . 

Although we find that neutron-rich winds from NDAFs 
are unlikely, several caveats should be discussed. Pruet et 
al. (2004) suggest that "bubbles" of neutron-rich material may 
escape the disk via chaotically-heated buoyant magnetic fila- 
ments, a picture similar to some models for GRBs (Narayan 
et al. 1992; Katz 1997; Kluzniak & Ruderman 1998). Al- 
though this possibility cannot be ruled out, current thin disk 
simulations do not find significant energy deposition in low 
density, coronal regions (e.g., Hirose, Rrolik, & Stone 2006), 
and whether such low-Y e bubbles can remain neutron-rich de- 
spite the pair-capture deneutronization that accompanies such 
chaotic heating is unclear. If the Blandford-Znajek mecha- 
nism or v— v annihilation above the disk's rotation axis pow- 
ers the GRB outflow instead of a disk wind, the base of the 
GRB-producing jet may be effectively baryon-free because 
the field lines would then thread the black hole's event hori- 
zon instead the disk midplane; what ultimately sets the wind's 
baryon-loading in this case is unclear. Our calculations show 

15 As discussed in i]4.4l an exception may arise for neutrino-driven winds 
produced at very early times after core bounce following the accretion- 
induced collapse (AIC) of a white dwarf (Dessart et al. 2006); similarly, in 
the absence of an overlying, accreting stellar mantle, an early-time ejection 
of ~ 10~ 3 - 1O~'M0 of low-Y e material may result from the core bounce 
shock's "break-out" or the subsequent neutrino-heated shock revival (e.g., 
Hillebrandt, Wolff, & Nomoto 1984; Fryer et al. 1999a) 



that neutrino-heated disk winds will form a modest entropy 
"sheath" around such a baryon-free jet. If, however, the wind 
encasing the jet possesses a much higher entropy (e.g., Pruet 
et al. 2001), cross-field neutron diffusion is more effective 
(Levinson & Eichler 2003; McKinney 2005b) and may result 
in asymptotically neutron-rich polar outflow. Furthermore, 
even if the outflow's field lines don't thread the disk, it is 
in principle possible that some form of chaotic mass-loading 
may pollute the baryon-poor base of the jet with matter from 
the neutron-rich disk midplane. We note, however, that cur- 
rent simulations find very little matter entraining the jet from 
the disk, thus requiring implementation of a numerical den- 
sity floor along the polar axis (e.g., Proga & Begelman 2003; 
Proga et al. 2003). 

6.3. Thick Disks 

Although NDAFs exist over a range of accretion rates 
that are relevant to both long and short-duration GRBs, 
GRB-producing outflows can also be powered by the accre- 
tion of matter that is not efficiently neutrino-cooled. Such 
geometrically-thick, radiatively inefficient accretion flows 
(RIAFs) exist for both M D < M ign ("low-M D RIAFs"), for 
which the density is too low for efficient neutrino cooling near 
tfisco, and for M D > M gn ("high-M RIAFs"), for which mat- 
ter advects into the black hole faster than it cools (e.g., Di 
Matteo, Perna, & Narayan 2002). HighM D RIAFs, while pos- 
sibly relevant to "prompt" collapsars (MacFadyen & Woosley 
1999), are probably most relevant to compact object mergers 
given their smaller expected disk radii and shorter accretion 
timescales. Low-Mo RIAFs are relevant to both the late stages 
of collapsars and compact binary mergers. 

Like NDAFs, high-Mo RIAFs enter /3-equilibrium be- 
fore accreting and, although they are not as neutron-rich as 
NDAFs, they also typically have Yf <C 0.5 (e.g., Surman & 
McLaughlin 2004; Lee et al. 2005, hereafter L05). High-Mo 
RIAFs are confined to radii in the disk smaller than the "trap- 
ping" radius R u the point interior to which matter has insuf- 
ficient time to cool before accreting; R t exists outside R lsco 

for mass accretion rates greater thanM t w 9(2)a l J\MQ s" 1 for 
a = 0(0.95) and M = 3M (CB07). Although radiatively inef- 
ficient on the whole, disks with Mo > M still release a sub- 
stantial neutrino luminosity because M t is large and accretion 
is still efficiently cooled for radii larger than R t ; indeed, for 
steady-state accretion, Di Matteo, Perna, & Narayan (2002) 
find that as Mo increases beyond M t , L v saturates to a value 
~ 10 53 ergs s" 1 . This substantial neutrino luminosity will 
drive significant neutrino-heated mass-loss, thereby severely 
limiting the asymptotic Lorentz factor V of any outflow driven 
from the RIAF's surface. Using the simulations of L05 (their 
Fig. 6), we estimate that the neutrino sphere temperature of a 
high-Mo RIAF near R t is roughly T v ~ 3 MeV (correspond- 
ing to a mean neutrino energy (e„) ~ 10 MeV); thus, from 
equation ( fT9] l we estimate that the minimum, neutrino-driven 
mass-loss rate from the innermost radii of high-Mo RIAFs is 
M t h « 10~ 3 - 10~ 2 M Q s" 1 . For a wind driven from the surface 
of a high-Mo RIAF to reach Y > 100 would thus require an 
MHD luminosity E > 10 53 - 10 54 ergs s" 1 , which is compara- 
ble to the entire accretion power for Mo ~ M t ~ 1 - 10M Q s" 1 
and is substantially larger than the outflow powers typically 
inferred from short-duration GRBs (e.g., Bloom et al. 2005). 
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GRB-producing outflows, whether neutron-rich or not, are 
thus unlikely to originate from disk winds produced by high- 
Mo RIAFs. 

Unlike NDAFs and high-Mo RIAFs, low-M D RIAFs may 
not enter /3-equilibrium on an accretion timescale (eq. fl28l ); 
thus, the electron fraction in a low-Mo thick disk's midplane 
can remain approximately equal to that of the material used 
to form the disk initially. Because of the low neutrino lumi- 
nosity and positive Bernoulli parameters of low-Mo RIAFs, 
viscous heating likely dominates neutrino heating in outflows 
driven from low-Mo thick disks, thereby making deneutron- 
ization by neutrinos unlikely. Although nondegenerate pair- 
captures may drive Y e — > 0.5 depending on the precise viscous 
energy deposition profile in the wind, this possibility is less 
likely than for NDAF outflows because /3-equilibrium is al- 
ready slow in the disk midplane and because the accretion and 
outflow advection rates for thick disks are comparable. Fur- 
thermore, if acceleration from the midplane is enhanced due 
to magnetocentrifugal slinging, pair-capture deneutronization 
would be further suppressed. Thus, outflows driven from low- 
Mo RIAFs probably retain the electron fraction of the disk's 
original composition. As a result, outflows driven at late times 
from collapsars, which are fed from the progenitor star's rel- 
atively neutron-poor envelope (e.g., Y e ~ 0.5 for a purely He 
composition), probably do not possess a significant neutron 
excess. In this case, significant 56 Ni could be produced in the 
disk's outflow, powering a bright SN (MacFadyen & Woosley 
1999). 

Although disks formed from compact object mergers are 
usually initially neutron-rich, they are compact and nearly 
all of the matter goes through an NDAF phase in which 
/3-equilibrium modifies the initial composition of the disk; 
their composition at late times therefore depends on how they 
transition from an NDAF to a low-Mo thick disk. Depend- 
ing on a, the size of the disk, and the mass of the black 
hole, the composition of the thick disk immediately follow- 
ing the NDAF to low-Mo RIAF transition at t lgn can be either 
neutron-rich or proton-rich (see eq. |29l ). Although by no 
means assured, neutron-rich GRBs outflows may thus be pos- 
sible from thick disks in both collapsars immediately follow- 
ing ?j gn (especially if the variability imposed by the NDAF- 
thick disk transition actually causes the GRB; Giannios 2007) 
and in short-duration GRBs from compact object mergers. In 
either case, neutron-rich GRB outflows should be restricted to 
events with accretion power (and thus maximum GRB lumi- 
nosity) less than L^ign (eq. Q). 

6.4. Non-Relativistic Neutron-Rich Winds 
6.4.1. Optical Transients 

Relativistic, GRB-producing outflows from proto- 
magnetars and hyper-accreting disks are possible under 
specialized circumstances, but non-relativistic winds are 
also likely to be present, probably occur in a wider variety 
of progenitors, and probably carry more total energy. The 
composition of such non-relativistic winds may also lead 
to an observable signature, most directly via radiation from 
ejecta that is reheated by the decay of radioactive elements 
in the wind. One possibility in the case of magnetar birth is 
that the heavily mass-loaded wind that emerges at early times 
(carrying a total mass up to ~ 0.01 -0. 1M Q for a millisecond 



rotator; Dessart et al. 2007; see MTQ07, Fig. 10) could 
produce a bright SN-like event. However, by combining the 
mass-loss rate that accompanies a significant neutron excess 
(eq. ffT8l ) with the fiducial PNS cooling evolution L v (t ) given 
in equation ®, we find that the total PNS mass-loss capable 
of being processed into 56 Ni under NSE (which requires no 
significant neutron excess; Hartmann, Woosley, & El Eid 
1985) cannot exceed M™ x - M±(/> n dt ~ 10 _3 M Q , much 
too small to contribute appreciably to an optical light curve 
powered by 56 Ni, and subsequent 56 Co, decay (e.g., see 
Kulkarni 2005, Fig. 7). 16 

The modest optical luminosity associated with the decay 
of a mass ~ M™ x of Ni is consistent with the rather strin- 
gent upper limits on the SN component accompanying some 
short-duration GRBs (e.g., GRB050509B; Hjorth et al. 2005); 
this supports the viability of a model in which short-duration 
GRBs are powered by the rapid spin-down of a magnetized, 
rapidly rotating magnetar, formed following the accretion- 
induced collapse (AIC) of a white dwarf (Usov 1992) or re- 
sulting from the merger of a double neutron star binary (a 
"super-pulsar"; Rosswog et al. 2003). It is less clear how 
much of the substantial mass-loss driven from the accretion 
disk formed during the merger of two compact objects will be 
processed into 56 Ni. Quantifying this will require additional 
work. 

Even if little mass is ejected with Y" > 0.5 (capable of pro- 
ducing 56 Ni), both proto-magnetars and hyper- accreting disks 
could in principle produce detectable transients due to the 
presence of neutron-rich non-relativistic outflows. For exam- 
ple, a proto-magnetar that is born rotating sufficiently rapidly 
to produce late-time, relativistic neutron-rich matter must also 
eject a total mass M n > 0.1M Q in slower (v > 0.3 c) free neu- 
trons at earlier times. The detectability of such neutron-rich 
non-relativistic outflows is, however, uncertain. Any neutrons 
that ultimately remain free in the wind will /3-decay into pro- 
tons on a timescale Tp 900 s at a radius > 10 13 cm; thus, 
one observable manifestation of neutron-rich outflow may be 
"Macronovae" powered by the thermal energy released by 
this decay («6 x 10 49 (M n /0.1M Q ) ergs), which may be de- 
tectable on hour-day timescales following the birth of the cen- 
tral object (Li & Paczyriski 1998; Kulkarni 2005). 

6.4.2. r-process Nucleosynthesis 

Rather than leaving free neutrons, the decompression of 
slowly moving, modest entropy, and moderately neutron-rich 
(0.1 < Y" < 0.4) matter is more likely to produce anomalous 
neutron-rich isotopes (e.g., 62 Ni, 56 Zn, 68 Zn, 87 Rb, 88 Sr; Hart- 
mann, Woosley, & El Eid 1985), and in some cases may be 
capable of producing r-process elements (e.g., Freiburghaus 
et al. 1999). Thus, as a final application of our calculations, 
we briefly consider the possibility of r-process element syn- 
thesis in proto-magnetar and NDAF winds. Given their in- 
trinsic neutron-rich nature, outflows from neutron star for- 
mation have long been considered one of the most promis- 
ing r-process sites (Woosley & Hoffman 1992; Meyer et 

16 Enhanced Ni production is still possible in the core collapse context via 
early energization of a successful SN shock due to rapid spin-down (TCQ04). 
A proto-magnetar's spin-down luminosity is substantially enhanced over the 
vacuum-dipole rate at early times following the launch of the SN shock by the 
excess magnetic flux opened by neutrino-heated, centrifugally-driven mass- 
loss (B06). In this case, the enhanced Ni yield is produced by additional 
shock heating of the stellar progenitor, not directly in the wind itself. 
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al. 1992); however, the conditions necessary for a successful 
third-peak r-process have not been realized in detailed studies 
of non-rotating, non-magnetized PNS winds (QW96; Otsuki 
et al. 2000; Sumiyoshi et al. 2000; Wanajo et al. 2001; T01). 
In MTQ07 we studied the effects of magnetic fields and rota- 
tion on the r-process in PNS winds for P > P n and constant 
Y e ; although we did not find solutions with a successful third- 
peak r-process (based on the criteria of Hoffman, Woosley, 
& Qian 1997, hereafter HWQ97), we did not consider the 
additional benefits of low Y" caused by very rapid rotation 

(P<Pnl 

Even accounting for the possibility of low Y e discussed in 
this paper, we find that all of the NDAF and most of the proto- 
magnetar wind solutions we have considered still fail to meet 
the criteria for third-peak r-process of HWQ97, mainly be- 
cause the beneficial effects of low Y e are counteracted by the 
detrimental effects that the accompanying rapid advection has 
on the wind's asymptotic entropy S". 17 One possibility is 
that energy deposition by acoustic or MHD waves could raise 
the asymptotic entropy of low-y c winds (Suzuki & Nagataki 
2005; MTQ07). In addition, despite their extremely low en- 
tropy, some of the most rapidly rotating and highly magne- 
tized proto-magnetar solutions that we have calculated (e.g., 
the B v = 10 16 G and = 9000 Hz solution shown in Figure 
do eclipse the third-peak r-process threshold of HWQ97. 
This occurs because r-process synthesis of nuclei with mean 
mass A (typically ~ 195 for the third peak) by neutron cap- 
tures on seed nuclei with mean proton number Z (typically 
- 30-40) is possible for Y° <Z/Am 0. 15 -0.20, even in out- 
flows with vanishingly small S a (see the discussion in Wheeler 
et al. 1998). The dynamical timescales for our very low Y" so- 
lutions are, however, much shorter than those considered by 
HWQ97 and several commonly made assumptions become 
suspect for such rapidly expanding outflows (Meyer 2002); 
thus, detailed nucleosynthesis calculations, which include all 
of the relevant rate equations, need to be completed in the 
regime of low S", low Y", and very rapid outflow before r- 
process success is assured. 

The very rapid rotation and strong magnetic fields asso- 
ciated with the successful r-process winds that we find are 
extreme and will not accompany most core-collapse SNe; 
the substantial amount of neutron-rich material (> 0.1M Q ) 
ejected by such a proto-magnetar, however, means that just a 
few very rapidly rotating proto-magnetar births of this kind 
per Myr would noticeably affect the Galactic r-process abun- 
dance (Qian 2000). Conversely, if the r-process yield from 
events like this do not resemble the observed abundances, 
the number of sub-millisecond magnetar births in our Galaxy 
can be strongly constrained; similar constraints can be placed 
on the incidence of the accretion-induced collapse (AIC) 
of a white dwarf, which produces a similar yield of low Y e 
material (e.g., Woosley & Baron 1992; Fryer et al. 1999a; 
Dessart et al. 2006, 2007). 

17 By the arguments given in [|2] for a PNS wind to remain neutron-rich 
the thermal energy deposited in the wind per baryon cannot exceed the neu- 
trinosphere's mean neutrino energy Thus, because most heating occurs 
near the PNS surface (at a temperature 0.7T„; QW96, eq. [47]), the en- 
tropy added to a neutron-rich wind cannot exceed AS ~ (e^/OJTjj Kb 5 
kg baryon - ' (compare S" and Y" in Tables [2] and [3J, which is substan- 
tially less than that deposited in non-magnetized, non-rotating PNS winds 
AS w GMmn/Q.lTvRv Ri 70(L o<!l 5i /8)- 1 / 4 k B baryorT 1 (QW96; T01). 
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TABLE 1 

Definitions for Commonly Used Variables 



Variable Definition 



Y e Electron fraction (ratio of free protons to nucleons) 

F e D Disk midplane electron fraction 

F c ° Electron fraction at the base of the wind 

Y° Asymptotic electron fraction in the wind 

Y^(r) Electron fraction in local weak equilibrium, defined via dY e /dt\ ^ =0(eq. Q4J) 

'e 

Y^ Asymptotic electron fraction in neutrino absorption equilibrium (Y" l (r — * oo) = Y"; eq. (6)) 

ya,m Asymptotic electron fraction that obtains in the co-rotating, strong B limit in our simulations (eq. 1 16| ) 

M Wind mass-loss rate 

Afth Purely thermal, neutrino-driven mass- loss rate (eqs. [19], (24l ) 

Mo Disk mass accretion rate 

Mign Minimum "ignition" NDAF accretion rate (eq. fTl ) 

Mf [ Total neutron mass diffusion rate into the accretion disk's polar region from an encasing baryon-rich wind (eq. [271 ). 

r Distance along the outflow to the center of the PNS or black hole 

R v Radius of the PNS neutrino sphere and the base of the PNS wind 

Distance from the black hole to the base of the NDAF wind (see Fig. [6} 

s Distance along the outflow to the monopole center in our NDAF wind calculations (see Fig.[6j 

so Distance from the monopole center to the base of the NDAF wind (see Fig. [61 

R s Black hole's gravitational radius (GM/c~, where M is the black hole mass) 

Risco Radius of the black hole's innermost stable circular orbit 

R\ S n 'Ignition" radius interior to which accretion proceeds through a thin NDAF instead of a thick disk (eq. 1 1 1 1 > 

R p Radius of the NDAF's peak integrated neutrino emission 

Rp Radius interior to which a thick disk enters (3— equilibrium on an accretion timescale (eq. |28| ) 

R n Radius interior to which a thick disk favors a neutron-rich composition in -equilibrium (eq. [291 ) 

O Rotation rate of the PNS and the base of the PNS wind 

On Rotation rate above which the PNS wind is significantly neutron-rich (Yf ^ 0.25; see eq. |16| ) 

£1% Keplerian rotation rate of the accretion disk 

a Magnetization (potential asymptotic Lorentz factor) of PNS/NDAF winds (eqs. 1131 , 1151 , 1211 ) 

d"mux Maximum magnetization of neutrino-heated NDAF winds (eq. [25) ) 

aJJ^ Maximum magnetization of neutrino-heated NDAF winds with a neutron excess (eq. |26| ) 
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TABLE 2 
PNS Wind Properties 







£2 


P 






y" 


M 


a 


5» ) 


10 51 ergs s' 1 


G 


s -1 


ms 


gcrrT 3 






MqS" 1 




kB baryorT 1 


8 


10 16 


1000 


6.3 


4 X 10 12 


0.03 


0.48 


1.9 x 10~ 4 


10 


55 




lO 16 


4000 


1.6 


4 X 10 12 


0.03 


0.47 


7 x ur 4 


40 


28 




10 1<s 


6000 


1.0 


3 x 10 12 


0.03 


0.41 


4 x ur 3 


17 


16 




10 


8000 


0.8 


1.7 x 10 12 


0.04 


0.26 


6 x ur 2 


1.9 


10 




10 


9000 


0.7 


1.3 x 10 12 


0.05 


0.19 


2.4 x 10~' 


0.6 


8 




10 15 


4000 


1.6 


4 x 10 12 


0.03 


0.47 


7 x 10" 4 


0.4 


28 




10 15 


6000 


1.0 


3 x 10 12 


0.03 


0.42 


4 x ur 3 


0.18 


17 




10 15 


8000 


0.8 


1.7 x 10 12 


0.04 


0.32 


4 x ur 2 


0.03 


1 1 




10 


9000 


0.7 


1.3 x 10 12 


0.05 


0.28 


9 x ur 2 


0.016 


10 




10 14 


4000 


1.6 


4 x 10 12 


0.03 


0.48 


7 x 10"* 


4 x 10~ 3 


35 




10 14 


6000 


1.0 


3 x 10 12 


0.03 


0.46 


1.8 x 10~ 3 


4 x ur 3 


23 




10 14 


8000 


0.8 


1.7 x 10 12 


0.04 


0.44 


7 x ur 3 


i.6 x ur 3 


18 




10 


9000 


0.7 


1.3 x 10 12 


0.05 


0.42 


1.5 x 10~ 2 


i.o x ur 3 


16 


3.5 


lO 16 


1000 


6.3 


7 x 10 12 


0.011 


0.50 


2.8 x 10~ 5 


70 


64 




lO 16 


4000 


1.6 


7 x 10 12 


0.011 


0.50 


1.0 x 10~ 4 


300 


31 




10 16 


6000 


1.0 


6 x 10 12 


0.011 


0.44 


6 x ur 4 


120 


17 




10 


8000 


0.8 


4 x 10 


0.014 


0.22 


0.016 


7 


9 




10 


9000 


0.7 


3 x 10 12 


0.015 


0.13 


0.10 


1.5 


7 




10 5 


4000 


1.6 


7 x 10 


0.011 


0.50 


1.0 x 10~ 4 


3 


31 




10 15 


6000 


1 


6 x 10 12 


011 


44 


6 x 10 -4 


1 2 


17 




10 15 


8000 


0.8 


4 x 10 12 


0.014 


0.24 


0.014 


0.08 


10 




10 15 


9000 


0.7 


3 x 10 12 


0.016 


0.18 


0.06 


0.024 


8 




10 14 


4000 


1.6 


7 x 10 12 


0.011 


0.50 


9 x 10~ 5 


3 x 10~ 2 


32 




10 14 


6000 


1.0 


6 x 10 12 


0.011 


0.46 


4 x ur 4 


1.5 x ur 2 


19 




10 14 


8000 


0.8 


4 x 10 12 


0.014 


0.37 


5xur 3 


2.6 x ur 3 


13 




10 14 


9000 


0.7 


3 x 10 12 


0.016 


0.34 


0.013 


1.2 x 10~ 3 


12 


1 


10 15 


7000 


0.9 


2.0 x 10 13 


0.003 


0.42 


9 x 10~ 5 


10 


14 




10 15 


8000 


0.8 


1.7 x 10 13 


0.003 


0.23 


1.0 x 10~ 3 


1.1 


9 




10 15 


9000 


0.7 


1.2 x 10 13 


0.004 


0.11 


0.011 


0.14 


6 



"Density at the base of the wind, set to enforce neutrino optical depth 7V ~ 2/3. 
*The asymptotic entropy of the wind. 
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TABLE 3 
NDAF Wind Properties 
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0.09 


0.32 


2.7 


0.016 


40 


10 




10 


60 


i j v , i a15 

1 .4 X 1 


2.2 X 10 


0.06 


0.47 


0.54 


8x 10 -4 


4 


1 3 




10 


60 


9x 10 14 


3 X 10 10 


0.04 


0.50 


0.09 


5 x 10 5 


0.5 


18 




LO 


70 


l.oX 10 


1.4X 10 


0. 1 1 


0.35 


1 .5 


0.028 


40 


1 1 




10 


70 


1.1 x 10 13 


2.0 X 10 


0.08 


0.46 


0.42 


l.Ox 10~ 3 


3 


13 




10 


70 


8 x 10 


3 X 10 


0.06 


0.50 


0.08 


6x 10~ 5 


0.5 


18 




10 


80 


5 X 10 


3 X 10 


0.05 


0.44 


0. 1 1 


0.39 


7 


10 




10 


80 


2.8 X 10 


-j v in 10 
i X 10 


0.05 


0.47 


0.08 


6x 10 3 


0.7 


1 8 




10 


80 


2.8 X 10 


4x 10 


0.04 


0.49 


0.07 


1.3 x 10~ 4 


0.2 1 
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10 15 


85 


1 9 X 10 14 


4 x 10 10 


04 


50 
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4 


ri 
J 


10 




10 M 


85 


1.1 x 10 14 


4x 10 10 


0.03 


0.50 


0.01 1 


0.04 


0. 18 


i 8 




10 13 


85 


6x 10 13 


4x 10 10 


0.03 


0.50 


0.009 


5x 10~ 4 


0.06 


24 




10 15 


87 


1 2x 10 14 


5 x 10 10 


03 


51 


5x1 3 


g 


3 


16 




10 M 


87 


8 x 10 13 


5 x 10 10 


03 


51 


5 x 10~ 3 


08 


12 


23 




10 13 


87 


2 4x 10 13 


6x 10 10 


03 


5 1 


5x1 3 


9 x 1 0~ 4 


029 


1L) 




in 15 
J u 


89 


1 .3 x 10" 


1 .U X 1U 


0.02 


0.5 1 


1 .5 x 10 3 


28 


3 


40 




1U 


vo. 

oV 


J X 1U 


1 .0 x 10' ' 


u.uz 


U.J 1 


l.OX LU 


u.z / 


n en 

\)A) I 






10 13 


89 


9xl0 12 


l.OxlO 11 


0.02 


0.51 


1.4 X 10" 3 


3xur 3 


0.008 


53 


0.01 


10 14 


50 


8xI0 14 


5xl0 9 


0.06 


0.11 


0.21 


2.1xl0~ 3 


40 


9 




10" 


60 


8xl0 14 


5xl0 9 


0.06 


0.12 


0.19 


2.2 xlO~ 3 


30 


9 




10" 


70 


5xl0 14 


5xl0 9 


0.06 


0.15 


0.10 


4xl0~ 3 


21 


9 




10" 


75 


1.3X10 14 


5xl0 9 


0.06 


0.27 


0.010 


0.04 


4 


12 




10" 


80 


4xl0 13 


5xl0 9 


0.06 


0.46 


6X10" 4 


0.7 


1.1 


17 




10" 


85 


2.2xl0 13 


5xl0 9 


0.06 


0.51 


8xl0" 5 


6 


0.6 


36 




10" 


87 


1.5 x 10 13 


5xl0 9 


0.06 


0.51 


2.6 x 10" 5 


17 


0.6 


55 




10" 


89 


1.3X10 13 


5xl0 9 


0.06 


0.51 


1.7 x 10" 5 


26 


0.6 


60 



"Azimuthal magnetic field at the base of the wind. 
''Density at the base of the wind. 

c The ratio of angular momentum lost in the wind to that lost through the disk (eq. 1221 ); solutions with Jw > Jd are unphysical. 
rf The asymptotic entropy of the wind. 
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TABLE 4 

Neutron Content of Outflows from GRB Central Engines 



Central hngine 


Neutron-Rich 
GRB Outflow? 


Conditions/Comments 


Magnetar, CC <a) 


Sometimes 


Subset with SN + long-duration GRB energy > 4 x 10 52 ergs; 
restricted to GRB luminosities < 3 x 10 49 ergs s" 1 ; see i| 4. 2 land ^6.11 

~ 0. IMq of non-relativistic (v > 0.3 c) free neutrons ejected prior to neutron-rich GRB outflow 


Magnetar, AIC (bl 


Sometimes 


Same as CC; 

additional ~ 10~ 3 — 10 _1 Mq non-relativistic low-7 e matter ejected at early times; ,c> 
SN-like component optically-dim due to < 10~ 3 Mq total 56 Ni production; see 36.4.11 


NDAF, CC 


Unlikely 


Disk midplane enters /3-equilibrium;' d) 

outflow enters neutrino absorption equilibrium (Y^ ~ Y^)\ see £|5.3| 


NDAF, COM (c) 


Unlikely 


Same as CC; 

additional ~ 10~ - 10 _1 Mq non-relativistic low-7 e matter may be ejected at early times' 1 ^ 


Low-M D RIAF (gl , CC 


Unlikely 


Disk midplane may not enter /3-equilibrium, and stellar mantle feeding the disk has Y e ~ 0.5; 
outflow likely viscously-driven with Ff ~ Y®\ see i|5.5|and *|6.3| 


I nw-Mr, RTAF COM 


Possible 


Neution star tidal debus feeding the disk has Y^ 0.5, but composition altered by /3-equilibnum at high A$u« 
low M/j-RIAF may remain neutron-rich during NDAF to RIAF transition; see fl5.5| 
due to low-Mo, probably accompanies only relatively long short-duration GRBs 


High-M D RIAF, CC 


Unlikely 


Disk midplane enters /3-equilibrium; 

v e -dominated radiation field accompanying deleptonization favors Y]f > 0.5; (h ' 
large neutrino-driven mass-loss likely precludes relativistic disk winds 


High-M D RIAF, COM 


Unlikely 


Same as CC 



"Core Collapse (CC) 
''Accretion-Induced Collapse (AIC). 

'e.g., Hillebrandt, Wolff, & Nomoto (1984); Fryer et al. (1999a); Dessart et al. (2006, 2007). 

''Pruet et al. (2003); B03a. 

'Compact Object Merger (COM). 

'e.g., Rosswog et al. (1999b). 

^Radiatively Inefficient Accretion Flow (RIAF). 

''e.g., Rosswog & Liebendorfer (2003) 



